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UNIVERSALITY IN ONE DIMENSIONAL HIERARCHICAL COALESCENCE 

PROCESSES 

A. FAGGIONATO, F. MARTINELLI, C. ROBERTO, AND C. TONINELLI 

Abstract. Motivated by several models introduced in the physical literature to study 
the non-equilibrium coarsening dynamics of one-dimensional systems, we consider a 
large class of "hierarchical coalescence processes" (HCP). An HCP consists of an infi- 

i^j nite sequence of coalescence processes {£'"'(•) } n >i- each process occurs in a different 

"epoch" (indexed by n) and evolves for an infinite time, while the evolution in subse- 
quent epochs are linked in such a way that the initial distribution of f ( n+1 ) coincides 
with the final distribution off'™'. Inside each epoch the process, described by a suitable 
simple point process representing the boundaries between adjacent intervals (domains), 
evolves as follows. Only intervals whose length belongs to a certain epoch-dependent 
finite range are active i.e. they can incorporate their left or right neighboring interval 

P^H with quite general rates. Inactive intervals cannot incorporate their neighbours and can 

Qh ' increase their length only if they are incorporated by active neighbours. The activity 

ranges are such that after a merging step the newly produced interval always becomes 
inactive for that epoch but active for some future epoch. 

Without making any mean-field assumption we show that: (i) if the initial distri- 
bution describes a renewal process then such a property is preserved at all later times 
and all future epochs; (ii) the distribution of certain rescaled variables, e.g. the domain 
length, has a well defined and universal limiting behavior as n — > oo independent of 
the details of the process (merging rates, activity ranges,...). This last result explains 
the universality in the limiting behavior of several very different physical systems (e.g. 

r^' ' the East model of glassy dynamics or the Past-all-model) which was observed in several 

simulations and analyzed in many physical papers. The main idea to obtain the asymp- 
totic result is to first write down a recursive set of non linear identities for the Laplace 
transforms of the relevant quantities on different epochs and then to solve it by means 
of a transformation which in some sense linearizes the system. 

p^ ' Keywords: coalescence process, simple point process, renewal process, universality, non- 

equilibrium dynamics. 
AMS 2010 Subject Classification: 60G55, 60B10, 82C24. 
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1. Introduction 

There are several situations arising in one dimensional physics in which the non- 
equilibrium evolution of the system is dominated by the coalescence of certain domains 
or droplets characterizing the experiment (e.g. large vapor droplets in breath figures 
or ordered domains in Ising and Potts models at zero temperature) which leads to 
interesting coarsening phenomena. As pointed out in the physics literature a common 
feature of these phenomena is the appearance of a scale-invariant morphology for large 
times. Many models, even very simple ones, have been proposed in order to capture 
and explain such a behavior (see e.g. llDBGH . IIDGY1I1 . IIDGY2I1 and ED). Supported by 
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computer simulations and under the key assumption of a well defined limiting behavior 
under suitable rescaling, physicists have derived some non trivial limiting distributions 
for the relevant quantities. 

In many cases the coalescence process dominating the time evolution has a hierar- 
chical structure which can, informally, be described as follows. 

Assume for simplicity that the state of the system is described by an infinite se- 
quence of adjacent intervals ("domains" in the physical language) with varying length 
and that its time evolution is governed by the merging of two consecutive intervals. 
Then there exist infinitely many epochs and in the n th -epoch only those domains whose 
length belongs to a suitable epoch-dependent characteristic range are active (or, bet- 
ter, n-active), i.e. they can incorporate their left or right neighbor interval with certain 
(bounded) rates which could depend on the epoch and on the length of the domain. 
Each epoch lasts a very long (mathematically infinite) time so that at the end of the 
epoch there are no longer n-active intervals provided that the total merging rate is 
strictly positive for any n-active domain. Then the next epoch takes over and the pro- 
cess is repeated. Clearly, in order for the successive coalescences to be able to eliminate 
domains created by previous epochs and therefore to increase the domain length, some 
assumptions about the active ranges should be made. If the n-th active range is the 
interval [d ffi ,, d^L) then we require that d^L = d^^. 

An interesting and highly non trivial example of a hierarchical coalescence process 
(HCP in the sequel) is represented by the high density (or low temperature) non- 
equilibrium dynamics of the East model after a deep quench from a normal density 
state (see llEJl ISBH for physical motivations and discussions and HFMRT1I1 for a math- 
ematical analysis) . The East model is a well known example of kinetically constrained 
stochastic particle system with site exclusion which evolves according to a Glauber dy- 
namics submitted to the following constraint: the 0/1 occupancy variable at a given 
site x£Z can change only if the site x + 1 is empty. In this case, if a domain represents 
a maximal sequence of consecutive occupied sites and if the particle density is very 
high, then the characteristic range of the length of active domains for the n <fc -epoch 
is [2 n_1 ,2 n ) and active domains can only merge with their left neighbor. Notice that 
with this choice for the active range the merging of two n-active domains automatically 
produces a n-inactive domain. This is a technical feature that will always be supposed 
true throughout the paper. 

Another interesting HCP is given by the Paste-all-model HDGY2H which was devised 
to model breath figures formed by coalescing droplets in one dimension. In this case all 
the domains are sub-intervals of the integer lattice, the n-active length interval is {n} 
and active domains merge with their left/right neighbor with rate one. 

In USED the authors, under the assumption that the scaled domain length has a well 
defined limiting behavior as n — > oo, computed the exact form of the limiting distribu- 
tion for the above defined HCP corresponding to East (see section C of USEH ). Under 
a finite mean hypothesis they find that the limiting behavior is exactly the same as the 
one computed in HDGY2H (always assuming the limiting behavior and the mean filed 
hypothesis) for the Paste-all-model, a fact that they describe as "surprising". 

Our main result, stated in Theorem 12.191 solves completely this enigma. In fact, 
without making any mean field hypothesis, we: 



ID HIERARCHICAL COALESCENCE PROCESSES 



(a) prove the existence of a well defined limiting behavior which is independent of the 
various merging rates; 

(b) classify the limiting distribution according to the initial one (i.e. the distribution at 
the beginning of the first epoch). 

Slightly more precisely the main content of our contribution can be formulated as fol- 
lows. Let £ denotes the random set of the separation points between the domains (do- 
main walls in physical jargon). Then, under very general assumptions on the merging 
rates and on the active ranges but always assuming d^lx = d^ hl for each n : 
(i) if at the beginning of the first epoch £ is described by a renewal point process (as 
implicitly done in the physical papers), then the same property holds for all times 
and all epochs; 
(ii) if Z^> denotes the domain length at the beginning of the n th -epoch rescaled by a 
factor l/rf ^ and if g^(-) denotes its Laplace transform, then g( n ~) — > g^ where 
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/-oo e -sx 

(s) = 1 — exp < — Co / dx > (1.1) 

provided that lim^o — s -^g^(s)/(l - g^\s)) = c (necessarily c G [0, 1]). More- 
over the above limit exists with cq = 1 when starting with a stationary renewal 
point process (which has therefore a finite mean) . If instead the initial law is in 
the domain of attraction of an a-stable law with a € (0, 1) then the limit exists 
with co = a. 
The above results, which can be generalized to exchangeable point processes, explain 
clearly why apparently very different physical systems (i.e. with different merging rates 
and/or active ranges) show the same asymptotic behavior. 

We want to stress here the crucial ideas behind the proof of our limit theorem. The 
first step goes as follows. Inspired by the form of the limiting distribution found by the 
physicists one uses the theory of complete monotone functions and Laplace transform, 
to show that for each n there exists a nonnegative Radon measure t^ on (0, oo) such 
that the Laplace transform for the n th -epoch, g^ n >, can be written as 

gW)(s) = 1 - exp ( - f e —t {n) (dx)\ . (1.2) 

1 7[l,oo) X ) 

Then one observes that the Laplace transforms {g^} n >i must satisfy a non-linear and 
highly non trivial recursive system of identities which, thanks to step one, translate into 
recursive identities for the measures t^ n >. In turn the latter can be solved to express the 
measure t^ in terms of £W in a simple form. Finally, the explicit form of *w allows us 
to pass to the limit n — >■ oo in the recursive identities and prove the main result. 

Coalescence processes (also called coagulation or aggregation processes) and their 
time-reversed analog given by fragmentation processes have also been recently much 
studied in the mathematical literature with different motivations and from different 
points of view (see e.g JA|], llBell and references therein). Most of the mathematical 
research focused on models with a certain mean-field character (i. e. the spatial position 
of the coalescing objects does not play any role) with some exceptions (see e.g HBe2ID . 
Although our model shows indeed a mean-field nature (see e.g. Remark [2 .161) due to 
the fact that the domain wall process £ is renewal or exchangeable at any future time t 
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if it was so at time t = 0, we have been able to explore some dynamical aspects of the 
HCP for which the geometrical alignment of the domains is relevant (see Section [3]). 

We conclude by mentioning that in HFMRT2H the methods developed here have been 
successfully applied to other HCP's, where a domain can also coalesce with both its 
neighboring domains as in HBDGH . In this class a particular interesting case is rep- 
resented by the model in which (roughly) the smallest interval merges with its two 
neighbors. In the mean field approximation and by forgetting how much physical time 
elapses between and during the merging events, one can derive a time evolution equa- 
tion for the domain size distribution in which the time variable t is a continuous ap- 
proximation of the discrete label n of the epochs. This equation has been rigorously 
analyzed in HGMH (see also llPell for an interesting review) by means of non-linear anal- 
ysis techniques. 
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2. Model and results 



In this section we introduce the main objects of our analysis, namely the simple point 
processes, the one-epoch coalescence processes and the hierarchical coalescence processes. 
Then we expose our main results. We start by recalling some basic notions of simple 
point processes, referring to IIDVI1 and HFKASH for a detailed treatment. 

2.1. Simple point processes. We denote by Af the family of locally finite subsets 
£ C E. Af is a measurable space endowed with the er-algebra of measurable subsets 
generated by 

{ZeAf : \£r\A 1 \=n 1 ,--- ,\£nA k \=n k }, 

Ai, . . . , Af. being bounded Borel sets in R and m,...,rik G N. On Af one can de- 
fine a metric such that the above measurable subsets correspond to the Borel sets 
HMKKH . We call domains the intervals [x,x'] between nearest-neighbor points x,x' 
in £ U {— oo, +00}. Note that the existence of the domain [—00, x'\ corresponds to the 
fact that £ is bounded from the left and its leftmost point is given by x' . A similar con- 
sideration holds for [x, 00]. Points of £ are also called domain separation points. Given 
a point x G M, we define 

4 := inf{£ >0: x-tG?}, d r x := inf{t > : x + t G £} , 

with the convention that the infimum of the empty set is 00. Note that if x G £ then d x 
(d x ) is simply the length of the domain to the left (right) of x. 

We recall that a simple point process (shortly, SPP) is any measurable map from 
a probability space to the measurable space Af. With a slight abuse of notation we 
will denote the realization of a SPP by £ while we will usually denote by Q its law on 
the measurable space Af. In what follows N (N+) will denote the set of nonnegative 
(positive) integers. 

Definition 2.1. 

(i) We say that a SPP £ is left-bounded if it has a leftmost point and has infinite cardi- 
nality. 

(ii) We say that a SPP £ is Z-stationary if £ c Z and its law Q is invariant by Z- 
translations, i.e. if for any x ETLthe random set £ — x has law Q. 

(Hi) We say that a SPP £ is stationary if its law Q is invariant under R-translations, i.e. 
if for any x G K the random set £ — x has law Q. 

Thanks to Theorem 1.2.2 in HFKASH and its adaptation to the lattice case, if £ is Z- 
stationary or stationary, then a.s. the following dichotomy holds: £ is unbounded from 
the left and from the right or £ is empty. In the sequel we will always assume the first 
alternative to hold a.s. and we will write £ = {x^ : k € Z} with the rules: xq ^ < xi 
and Xk < Xk+i for all ieZ. In the case of a left-bounded SPP, we enumerate the 
points of £ as {x/c; k G N} in increasing order. 

Remark 2.2. If ^ is Z-stationary and a.s. nonempty, then Q(0 G £) > and therefore the 
conditional probability Q(-\0 G £) is well defined. On the other hand, if £ is stationary, 
then Q(0 G £) = 0, the above conditional probability is therefore not well defined and 
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has to be replaced by the Palm distribution associated to Q IIDVI1 . HFKASH . We recall that, 
given the law Qofa stationary SPP with finite intensity 

A S :=E s (|en[0,l]|) 

and such that £ is nonempty Q-a.s., the Palm distribution Qq associated to Q is defined 
as the probability measure on the measurable space M such that 

Qo(A) = (1/Aq)E q ( \{x G f n [0, 1] : t x £, g A}\ ) , VAcN measurable 

(see Section 1.2.1 in HFKASH ). Trivially, Qo has support in 

A/- °° := {£ e A/" : og£, |en(-oo,o]| = |^n[o,oo)| = oo}. (2.1) 

Moreover Qo univocally determines the law Q since it holds that 

E G [/(e)] = A Q E Qo / m-t)dt 

for any nonnegative measurable function f onM (cf. Theorem 1.2.9 in HFKASL Theorem 
12.3.11 in IIDVII ). Notice that, by taking f = 1, one gets Xq = 1/Eq (xi). Consider now 
the space (0, oo) z endowed with the product topology with Borel measurable sets. Setting 
d k (£) = x k (Q - Xfc_i(0 for k G Z and £ G M^, the map A/" °° 9(^ (0,oo) z !5 a 
measurable injection, with measurable image. In particular, the Palm distribution can be 
thought of as a probability measure on (0, oo) z . As stated in Theorem 1.3.1 in HFKASH . a 
probability measure Q on (0, oo) z is the Palm distribution associated to a stationary SPP 
with finite intensity and a.s. nonempty configurations if and only if Q is shift invariant 
and its marginal distributions have finite mean. 

We now describe the main classes of SPP's we are interested in. 

Definition 2.3. Given a probability measure p, on (0, oo), we say that £ is a renewal SPP 
containing the origin and with interval law p, and write Q = Ren(p \ 0), if 

(i) G £; 

(ii) £ is unbounded from the left and from the right and, labeling the points in increasing 
order with xq = 0, the random variables d k = x k — Xfe-i, k G Z, are Ltd. with 
common law p. 

Definition 2.4. Given probability measures v and p, on M. and (0, oo) respectively, we 
say that £ is a right renewal SPP with first point law v and interval law p, and write 
Q = Ren(u, p), if 

(i) £ = { Xk , k G N} is a left-bounded SPP, 

(ii) the first point x$ has law v, 

(iii) dk = Xk — Xk-\ (k G N + j has law p, 

(iv) The random variables xq, {dk}ken + are independent. 

Definition 2.5. Given a probability measure p on N + with finite mean, we say that £is a 
Z-stationary renewal SPP with interval law p, and write Q = Renz(p), if 

(i) £ is Z-stationary and a.s. nonempty, 

(ii) w.r.t. the conditional probability Q(-|0 G £) the random variables dk = Xk — Xk-i, 
k G Z, are Ltd. with common law \x. 
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A basic example is the following. Consider a Bernoulli product measure on {0, 1} Z 
with parameter p. Any realization (Xj)j 6 z can be identified with the subset £ = {i e 
Z : Xi = 1}. The resulting SPP is a Z-stationary renewal SPP with geometric interval 
law. 

Remark 2.6. As proven in Appendix [C] given a probability measure p on N+, the law 
Q = Renz(fi) is well defined iff p has finite mean. Other properties of Z-stationary 
renewal SPP's are also discussed there. 

Definition 2.7. Given a probability measure p on (0, oo) with finite mean, we say that £ 
is a stationary renewal SPP with interval law p, shortly £ = Ren(p), if 

(i) £is a stationary SPP with finite intensity and £ is nonempty a.s., 
(ii) the random variables dk = Xk — x^-i, k G Z, are Ltd. with common law p w.r.t. 
the Palm distribution associated to Q. 

A classical example of stationary renewal SPP is given by the homogeneous Poisson 
point process, for which the interval law is an exponential. 

Remark 2.8. A stationary renewal SPP with interval law p having infinite mean cannot 
exists (see Proposition 4.2.1 in IIDVII ). As discussed after Theorem 1.3.4 in HFKASH . Q = 
Ren(p) if and only if the following holds: the random variables dk = Xk — Xk-i, k ^ 1, 
are Lid. with law p and are independent from the random vector (xq, x\), which satisfies 

POO 

Q(-x > u, xi > v) = X Q I (1-F(t))dt, F(t) := p((0,t\) , u,v>0. (2.3) 

J u+v 

We conclude with the definition of two large classes of "exchangeable" point pro- 
cesses. 

Definition 2.9. We say that £is a left-bounded exchangeable SPP containing the origin 
if 
(i) £ = {xk, k € N} is a left-bounded SPP containing the origin; 

(ii) Q, thought of as probability measure on (0, oo) N + by the map £ — > (xk — Xk-i ■ k G 
N_|_), is exchangeable (i.e. invariante under permutations HDllKll ). 

Definition 2.10. We say that £ is a stationary exchangeable SPP if 

(i) £is a stationary SPP with finite intensity and £ is nonempty a.s.; 
(ii) the Palm distribution Qq, thought of as probability measure on (0, oo) z by the 
map £ — > (xk — Xk-i : k G Z), is exchangeable. 

Remark 2.11. Any left-bounded or stationary renewal SPP is also exchangeable. 

2.2. One-epoch coalescence process. We describe here the class of coalescence pro- 
cesses which will represent the modular unity of the, yet to be defined, hierarchical 
coalescence process (HCP). For a reason that will become clear in the next section, we 
call it one-epoch coalescence process (OCP). 

This process depends on two constants < d m ; n < d max and on nonnegative bounded 
functions A^, A r defined on [d m i n , oo] which, with \(d) := Xe(d) + X r (d), satisfy the fol- 
lowing assumptions: 

(Al) X(d) > if and only if d € [d min , d max ), 
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(A2) if d, d! > d m i n , then d + d! ^ <i m ax- 
Trivially, (A2) is equivalent to the bound 2d m ; n ^ d max . 

The one-epoch coalescence process is a Markov process with state space J\f(d m i n ) 
given by the configurations £ G J\f having only domains of length not smaller than 

A/^min) = {£ G AA : 4 > d min , < > d min Vx G £} • (2.4) 

The stochastic evolution is given by a jump dynamics with cadlag paths (£(£) : t ^ 0) in 
the Skohorod space D([0, oo), J\f(d m [ n )) (cf. HB]]) and at each jump a point is removed. 
Formally, the Markov generator of the coalescence process is given by 

£/(£) = E( A ^i) + MO) [/(£ \ W) - /«)] • (2.5) 

We will write Pq for the law onD([0, oo), A/"(d m m)) of the one-epoch coalescence pro- 
cess with initial law Q on J\f(d m { n ) and Q t for its marginal at time t. 

We keep the discussion of the Markov generator at a formal level, since we prefer 
to give a constructive definition of the coalescence process. Here we give two rough 
alternative descriptions of the dynamics as a random process of points or as a random 
process of domains (intervals). 

Point dynamics. To each point x G £(0) we associate two exponential random variables 
T x g and T X)T of parameter Xe(d x ) and X r (d x ), respectively. We stress that d x and d x 
refer to the configuration £(0): at time the domains on the left and on the right of 
the point x are respectively [x — d x ,x] and [x,x + d x ]. All random variables must be 
independent. If t = T x ^ ^ T x>r and at time t— the point x — d x still exists, then we 
set £(*) = £(t— ) \ {x}. Moreover, we say that the two domains having x as separation 
point merge or coalesce at time t, and that the domain on the left of x incorporates the 
domain on the right of a; at time t. If t = T x>r < T x ^, and at time t— the point x + d x ' still 
exists, then we set £(£) = £(£— ) \ {x}. Moreover, we say that the two domains having 
x as separation point merge or coalesce at time t, and that the domain on the right of 
x incorporates the domain on the left of x at time t. Finally, if t = T x ^ r A T Xi i but the 
above two cases do not take place, then we set £(£) =£(£—). See Figured] below for an 
illustration. 

In order to formalize the above construction, we proceed as follows. Given t > 0, 
we define T t as the set of points x G £(0) such that T x j A T xx ^ t. On the set T t we 
define a graph structure putting an edge between points x, y G T^ if and only if x and y 
are consecutive points in £(0). Since the functions A^, A r are bounded from above, a.s. 
for any fixed time t the above graph T t has only connected components (clusters) of 
finite cardinality. Then, £(0) \ T* is included in £(s) for all s G [0, t], while the evolution 
of (£(s) : s G [0,t]) restricted to each cluster of T t follows the rules stated at the 
beginning, which are now meaningful a.s. since clusters have finite cardinality a.s. 

Domain dynamics. We give here only a rough description of the dynamics. In Section 
13.11 we will discuss in detail a basic coupling leading to the definition on the same 
probability space of the domain dynamics for all initial configurations £(0) G A/"(d m ; n ). 
One assigns to each domain A = [x, x'\ with length d present in £(0) an exponen- 
tial random variable TX of parameter \(d) and a coin Ca with faces —1,1 appearing 
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with probability \ r (d)/X(d) and \^(d)/\(d), respectively. All random variables must be 
independent. If t = Ta and if at time £— the domain A still exists, then at time £ the 
domain A incorporates its left domain (i.e. £(£) = £(£— ) \ {x}) if Ca = — 1, while A 
incorporates its right domain (i.e. £(£) = £(£— ) \ {x'}) if Ca = 1. 

We can now explain the dynamical meaning of assumptions (Al) and (A2): 

• (Al) means that a domain is active, i.e. it can incorporate another domain, iff 
its length d lies in [d min , <i max ). 

• (A2) means that a domain resulting from a coalescence is not active. 

As consequence, the following blocking effect appears: given three nearest-neighbor 
inactive domains Ai, A2, A3, the intermediate domain A2 is frozen, in the sense that 
its extreme points cannot be erased, see Figure [TJ 
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Figure 1 . An example of the one-epoch coalescence process starting 
from £(0). At time £ = 0, the domain of length d is inactive since 
d ^ d max . At time t\, site x disappears and since £1 = T X) e, the domain 
on the left of x incorporates the domain on the right of x. Analogously at 
£2 point y disappears since the domain on the right of y incorporates the 
domain on the left of y. The domain Ai and A3 are inactive since they 
are resulting from a coalescence. The domain A2 is frozen for £ > £2, 
due to the presence of Ai and A3. This illustrate the blocking effect. 

By definition of the one-epoch coalescence process, points can only be removed. 
Therefore, on any finite interval /, £(£) n I converges as £ — > 00 and the following 
Lemma follows at once. 

Lemma 2.12. For any given initial condition £ G AA(d m ; n ) the following holds: 

(i) £(£)c£( S )i/ S ^£, 

(ii) the configuration £(£) is constant on bounded intervals eventually in t, 
(Hi) there exists a unique element £(00) in M{d raax ) such that £(£) (1 1 = £(00) n I for 
all large enough t (depending on I) and all bounded intervals I. 

Due to the above Lemma, £(00), the SPP representing the asymptotical state of the 
coalescence process, is well defined. Our first main result is given by the following 
two theorems. It states that, starting from a left-bounded renewal (respectively a Z- 
stationary or stationary renewal) SPP £, at a later time £ the coalescence process £(£) 
remains of the same type. Moreover, there exists a key identity between the Laplace 
transform of the interval law at time £ = and time £ = 00. This equation, that we 
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call one-epoch recursive equation, will play a crucial role in a recursive scheme for the 
hierarchical coalescence process. 

Theorem 2.13 (Renewal property). Let u, /i be probability measures on R and [d m [ n , oo), 
respectively. Then, for all t G [0, oo] there exist probability measures vt,^t on R and 
[d m ; n , oo) respectively such that vq = v, fiQ = fi and: 

(V If Q = Ren(v, p) then Q t = Ren(i/ t , fx t ); 

(ii) ifQ= Ren(fi) then Q t = Ren(pt); 
(Hi) if Q = Ren z (fi) then Q t = Ren z (fi t ); 
(iv) IfQ = Ren(5 ,fi) then Q t (- | G £) = Ren(5 , fit); 

(v) limt^oo v t = Voo and lim^oo fi t = /j,^ weakly. 

Theorem 2.14 (Recursive identities). Let u, /i be probability measures on R and [d m m , oo), 
respectively, and let v t ,fj, t be the probability measures introduced in Theorem \2. 1 3\ 
(i) Consider the Laplace/characteristic functions 



G t (s) = / e- sx nt(dx) , sel+Uil, (2.6) 

J Knin,°o) 



H t (s) = / e- sx n t (dx) , sel+Uil. (2.7) 

Then, for any s G R + U iR, the following one-epoch recursive equation holds: 

1-G? 00 (s)= [1 - G (s)]e Ho ^ . (2.8) 

(ii) Consider the Laplace/characteristic function 



L t (s) = / e~ sx v t (dx) , s G R+U G 



(a) If X r = 0, then v t = vofor all t ^ 0. Hence L t (s) = Lo(s)for all t ^ 0. 

(b) If Xe = j\ r for some 7 G [0, 00), then, for any s G R+ U iM., 

Loo(s) = L Q (s) exp <^ — } . (2.9) 

I 1+7 J 

Moreover, if Q = Ren(u, p) 

gfltg) 

P S (x (0)G^(oo)) = e-^+T , (2.10) 

where cco(0) denotes the first point of the initial configuration £(0). 

The proof of Theorem l2.13l and Theorem 12 .141 is given in Section [3] and Section @] 

Remark 2.15. In (ii) we have analyzed two cases ((a) and (b)) motivated by the East 
model and by the Paste-all model. The arguments used in the proof of point (ii) could 
however be applied to other cases as well. We stress that the Laplace transform L t (s), 
s G R+, could diverge since u t has support on R. Therefore, the above identities in point 
(ii) have to be thought of as identities in the extended space [0, 00]. 

We point out that the one-epoch recursive equation (I2.8D univocally determines Hoo 
when knowing /io, d m [ n , d max . In particular, these three elements are the unique traces 
of the dynamics that asymptotically survive. In other words the precise form of the 
rates A^ and A r is irrelevant. In the case of a left-bounded renewal point process the 
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limiting first point law v^ does not share such a universality although the trace of \i 
and A r on v^ is only partial. 

Remark 2.16. Assume for simplicity that \x is concentrated on N + , so that the domains 



have integer length at any time. After properly constructing the Markov generator (I2.5D 
one could prove that 

d-l 

d t nt{d) = -X(d)fi t (d) + Y, [Hx) + K(d- x)] nt(x)fi t (d - x) . (2.11) 

x=l 

Note that if d is active then only the first addendum in the r.h.s. is present, while if d 
is inactive this first addendum is absent. From this observation, one easily obtains that 
dtGt = (1 — Gt)dtH t , and therefore 

l-G t (s) = (l-G (s))exp{H (s)-H t (s)}, Vt,a>0. (2.12) 

Taking the limit t — > oo one gets (I2.8D . This strategy has been applied in IISEI1 . where the 
treatment is not rigorous, and will be formalized in IIFMRT2I1 in order to treat other coales- 
cence processes as in IIBDGII . It could be applied to derive fl2.9D . While the Smoluchosw ski- 
type equation (I2.11D has a mean-field structure (see e.g. OAJX in proving (I2.8D and Q2.9D 
we have followed here a more constructive strategy and we have investigated how a do- 
main of given length can emerge at the end of the epoch or how a given point can become 
the first point for the configuration £(oo) at the end of the epoch. 



e (3) (0)=e (2) (oo) 



epoch 2 

h — i — , — • — i — , — • — h- e (2) (o) = e (1) (oo) 



epoch 1 

h , — • • • — i — m—m- e (1) (o) 



Ai of length 1 A 2 of length 3 A 3 of length 2 

Figure 2. An example of HCP dynamics, with d^ 1 ' = n. The distances 
between the points are, from left to right, 1, 1 (corresponding to Ai), 
2, 1, 3 (corresponding to A 2 ), 1, 1, 2 (corresponding to A3)... At the 
beginning of epoch 1, only the domains of length in 1 are active. In 
particular, Ai is active while A 2 and A3 are inactive. At the end of 
epoch 1, there are no more domains of length less than 2 (see Lemma 
I2.12D . At the beginning of epoch 2, domains of length 2 are active and 
at the end, there are no more domains of length less than 3, and so on. 
Note that an inactive domain as A2 can increase its length. 

2.3. The hierarchical coalescence process. We can finally introduce the hierarchi- 
cal coalescence process (HCP) . The dynamics depends on the following parameters and 



12 A. FAGGIONATO, F. MARTINELLI, C. ROBERTO, AND C. TONINELLI 



functions: a strictly increasing sequence of positive numbers {d^} n >i and a family 

of uniformly bounded functions xf 1 ' , \\. : [d^ n \oo] — > [0,A],n ^ 1. Without loss of 

generality we may assume that d^ = 1. We set as before A( n ) := A^ + X r n and we 
assume 

(Al) for any n e N+, X^{d) > if and only if d € [d^ n \ d^ n+1 ^); 

(A2) for any n e N+, if d, d' > d^ n \ then d + d' ^ d^ n+1 ^ (i.e. 2d^ ^ d ( - n+1 ^; 

(A3) lim^oo d(») = oo. 

For example one could take d^ = n or d^> = a™ -1 with a £ (1,2]. 

The HCP is then given by a sequence of one-epoch coalescence processes, suitably 
linked. More precisely, the stochastic evolution of the HCP is described by the sequence 
of paths {^ n \-)} n >i where each £( n ) is the random path describing the evolution of 
the one-epoch coalescence process with rates A^ , X r n , active domain lengths d^( n = 

d {n \ d { ™L = d( n+1 ) and initial condition £( n )(0) = ^""^(oo), n > 2. Informally we 
refer to £("■) as describing the evolution in the n' h -epoch. See Figure [2] for a graphical 
illustration. 

Theorem 12.131 gives us information on the evolution and its asymptotic inside each 
epoch when the initial condition is a SPP of the renewal type. If e.g the initial dis- 
tribution Q for the first epoch is Ren(i>, fj) we can use Theorem 12.131 together with 
the link ^ n+1 \0) = £( n )(oo) between two consecutive epochs to recursively define the 



measures uy 1 ', z/ n ) by 



/* (B+1) =/&\ » {1) =» 



!/(«+!) =!>), U^=U (2.13) 



With this position it is then natural to ask if, in some suitable sense, the measures 
/jW,i/W have a well defined limiting behavior as n — > oo. The affirmative answer is 
contained in the following theorem, which is the core of the paper. Before stating it we 
need a result on the Laplace transform of probability measures on [1, oo]. 

Lemma 2.17. Let fibe a probability measure on [1, oo) and let g(s) be its Laplace trans- 
form, i.e. g(s) = J e~ sx fi(dx) , s > 0. 

i) If 

lim--^f- = c , (2.14) 

sio 1 — g(s) 

then necessarily < cq < 1; 
ii) The existence of the limit (I2.14D holds if: 

a) fj, has finite mean and then cq = 1 or 

b) for some a € (0, 1) p, belongs to the domain of attraction of an a-stable law or, 
more generally, /x((x,oo)) = x~ a L(x) where L(x) is a slowly varying 1 function at 
+oo, a € [0, 1], and in this case cq = a. 

Remark 2.18. One could wonder if the limit (I2.14D always exists. The answer is negative 
and an example is given in Section [B] 

A function L is said to be slowly varying at infinity, if for all c > 0, lim L(cx) / L(x) — 1. 
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The proof of Lemma [2 .171 is discussed in Appendix [Al 

Theorem 2.19. Let v, p, be probability measures on R and [1, oo) respectively and let g(s) 
be the Laplace transform of ' \x. Let Q be the initial law of^ and suppose that Q is either 
Q = Ren(u,p) or Q = Ren{p) or Q = Renz(fi). For any n € N + let X( n ' be a random 
variable with law ^ defined in (127131) so that g(s) := E[e" sX(1) ]. 

If (I2.14D holds for g then the rescaled variable Z^ n > := Xw /d( n ) weakly converges to 
the random variable Z^°°^ = Z^° whose Laplace transform is given by 

9^' 



(s) = 1 — exp < —cq dx>. (2.15) 

The corresponding probability density is of the form z Co (x)I x j> i where z CQ is the continu- 
ous function on [1, oo) given by 



oo 



= 1^ ii - Pk(x)l x>k , (2.16) 

fc=i 



where p±(x) = 1/x and 



Pk+i(x) = / dxi--- dx k —r TT _ 

J\ J\ x- }^ i=l Xi j=1 Xj 



k > 1 



The proof of Theorem 12.191 is discussed in Section [6721 

Remark 2.20. The remarkable fact of the above result is that the only reminiscence of 
the initial distribution in the limiting law is through the constant cq which, as proved in 
Lemma \2.14\ is "universal" for a large class of initial law /i. Hence the term universality 
in the title. We also stress that, starting with a stationary or Jj-stationary renewal SPP, 
the weak limit of Z^ always exists and it is universal (cq = 1), not depending on the 
rates \f\ \ { r n) . 

,(oo) 



Remark 2.21. We point out that the asymptotic Laplace distribution g K Co ' can be written 
also as 



'(s) = 1 — exp < — cq / dx > = 1 — exp < — coEi(s) > 



where Ei(-) denotes the exponential integral function. This is indeed the form appearing 
in HDGY2H and USEH with cq = 1 (see previous Remark). 

If the law fi has finite mean then by the above Theorem combined with ii) of Lemma 
|2.17l we know that Z^ weakly converges to the random variable ZJ . Actually we 
can improve our result to higher moments. 

Proposition 2.22. In the same setting of Theorem \2. 1 9\ assume that d^ = a n ~ x for 
some a € (1,2], and that fi has finite k-th moment, k e N+. Then, for any function 
f : [0, oo) — > M. such that \f(x)\ ^ C + Cx k for some constant C, it holds 

lim E[/(2< n ')] =E[f(z[ oo) )]. (2.17) 
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Remark 2.23. The choice d^ = a n_1 in Proposition \2.22\ is technical and could be 
relaxed, but at the price of extra hypotheses (that would not include the case d^> = n, for 
example). In order to keep the computations as simple as possible we decided to focus on 
this particular example which is of interest for applications to the East model. 

The proof of Proposition 12.221 can be found in Section 16.51 Next we concentrate 
on the asymptotic behavior of the first point law when starting with a left-bounded 
renewal SPP. 

Theorem 2.24. Let u, p be probability measures on E and [1, oo) respectively and consider 
the hierarchical coalescence process such that the initial law Q of^ is Ren(u, p). Assume 

Aj n) = 7 A(. n) , Vn>l, (2.18) 

for some 7 £ [0, 00), and let, for any n e N+, Xq be the position of the first point of the 
HCP at the beginning of the n-th epoch. 

If the limit (I2.14D exists for the Laplace trasform g of p then, as n — >■ 00, the rescaled 
random variable yW := Xq/S- 1 "^ weakly converges to the positive random variable 
Ycq with Laplace transform given by 

-sY™h _..„ f c f l-e~ s y 



E(e^ c o ) = exp i — / dy } , s e R+ (2.19) 

[ 1 + 7 7(0,1) V J 

We point out that if \\ n ' = for all n ^ 1, the first point does not move 2 . In 
particular, its asymptotic is trivial. Theorem |2.24l is proven in Section [631 
Finally, we evaluate the surviving probability of a given point: 

Theorem 2.25. Let u, p, be probability measures on K and [1, 00) respectively and consider 
the hierarchical coalescence process with initial law Q. Assume 

\^ = 7 A(") , Vn > 1 , (2.20) 

for some 7 € [0, 00) and let, for any n € N +J , Xq be the position of the first point of the 
HCP at the beginning of the n-th epoch. 

If the limit (I2.14D exists for the Laplace transform g of p then, as n — > 00, 

(i) if Q = Ren(v, pi), then 

p s (xJ m) = xi 1) ) = (i/dW)^ (1+o(1)) , 

(ii) if Q = Ren(p I 0), then 

p s (oe^(o)) = (i/^)) co(1+o(1)) 

Note that (ii) does not depend on the value of 7. Theorem [2^25] is proven in Section 

[Ml 

Extension of the above results to one-epoch coalescence process or hierarchical co- 
alescence process with initial law Q describing an exchangeable SPP will be discussed 
in appendix [0 



9 
This is the case for the HCP associated to the West version of the East model, i.e. to the kinetically 

constrained model with Glauber dynamics for which the occupation variable at x can be updated iff x — 1 

is empty 
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3. Renewal property in the OCP: Proof of Theorem 12.131 

In this section and in the next one we will prove our results concerning the one- 
epoch coalescence process (Theorem 12.131 and Theorem I2.14D in a more general set- 
ting, namely when the interval [d m i n ,d max ) is replaced by a more general set A C 
(0, oo). More precisely, let A^, A r be bounded nonnegative functions on (0, oo] and set 
A = Xe + X r . We assume that 

(Al') A(d) > if and only if d G A, 
(A20 if d, d! > inf(A), then d + d! <£ A. 

Above, d mm ■= inf (A) denotes the inflmum of the set A. When A = [d m in, d max ), (Al') 
and (A2') coincide with assumptions (Al) and (A2), respectively. A domain is called 
active if its length belongs to A. The initial distribution Q of the one-epoch coalescence 
must be supported in [inf (A), oo). In (I2.6D and (I2.7D the integration domains become 

[inf (.A), oo) and A, respectively. 

The proof of Theorem 12.131 requires the definition of a universal coupling, i.e. the 
construction on the same probability space of all one-epoch coalescence processes ob- 
tained by varying the initial configuration. This coupling will be relevant also in the 
proof of Theorem l2.25l (ii) . 

3.1. Universal coupling for the domain dynamics. In Section [2] we have introduced 
some enumerations of the points in £ G J\f, depending on the property of £ to be 
unbounded both from the left and from the right, or only from the left. It is convenient 
here to have a universal enumeration. To this aim, given £ G Af, we enumerate its 
points in increasing order, with the rule that the smallest positive one (if it exists) gets 
the label 1, while the largest non-positive one (if it exists) gets the label 0. We write 
N(x, £) for the integer number labeling the point x € £. This allows to enumerate the 
domains of £ as follows: a domain [x,x'] is said to be the k th -domain if (i) x is finite 
and N(x, £) = k, or (ii) x = — oo and N(x', £) = k + 1. Recall that if x = — oo, then £ is 
unbounded from the left and x' is the smallest number in £. 

We set ||A||oo = snp de ^\(d) an< ^ we define A^ = A r , A* = Xi. Obviously A = 
Xg + A r = X\ + A* . This change of notation should help the reader. Indeed, in the point 
dynamics a point x is erased by the action of its left (right) domain of length d with 
rate Xi{d) (X r (d)). On the other hand, as explained again below, if we formulate the 
model in terms of a domain dynamics then a domain of length d disappears because of 
the annihilation of its left (right) extreme with probability rate X\ (d) ( A* (d) ) . 

We consider now a probability space (fl, F, P) on which the following random ob- 
jects are defined and are all independent: the Poisson processes T^ = {Tm ' : m G N} 
and T^ = {Tm : m 6 N} of parameter HAHoo, indexed by k G Z, and the random 
variables Um , Um , uniformly distributed in [0, 1], indexed by k G Z and m G N. 

Next, given ( G M(d m [ n ) and oj G O, to each domain A that belongs to ( we associate 
the Poisson process T^ if A is the k-th domain in (. In this case, we write T^ A ' instead 

of T {k) . Similarly we define f (A) , U^\ U^ ] . We define W t [u, C] as the set of domains 
A in C such that 

Is G [0,t] : s G T (A) uf (A) , or s G r (A,) uf (A ' } for some domain A'neighboring a| ^ 
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On Wt[w, C] we define a graph structure putting an edge between domains A and A' if 
and only if they are neighboring in (. Since the function A is bounded from above, we 
deduce that the set 

£>(£) := \oj : Wi[w, C] has all connected components of finite cardinality Vt ^ > 

has P-probability equal to 1. Note that the event B(() depends on ( only through the 
infimum and the supremum of the set {N(x,() 6 2 : a; G (}. By a simple argument 
based on countability we conclude that P(B) = 1, where B is defined as the family 
of elements u> G U belonging to n^ e ^/-(^ min )0(C) and such that all the sets T^[uS\ and 
f (fc) M, k G Z, are disjoint. 

In order to define the path {£(s)} s >o = {^(s,w)} s >o we first fix a time t > and 
define the path up to time t. If uj <^B, then we set 

£(s) = C, Vs€[0,t]. 

If uj G B, recall the definition of the graph Wt[w, (]. Given a set of domains V we write 
V for the set of the associated extremes, i.e. x G V if and only if there exists a domain 
in V having x as left or right extreme. Moreover, we write Vt[w,C] for the set of all 
domains in ( that do not belong to Wt[w, (]. We define 



e(«)nV t [w,C]:=Vt[a;,C], VsG[0,t], (3.1) 



i.e. up to time t all points in Vt[co, C] survive. Let us now fix a cluster C in the graph 
W t [co, (]. The path (£(s) nC : s G [0, £]) is implicitly defined by the following rules (the 

definition is well posed since uj G B). If s G [0, i] equals T„ with A = [x, x'] G C and 
x, x' G £(s— ), then the ring at time T„ is called legal if 

t4 A) < ^if,,"^ (3.2) 



and in this case we set£(s)nC := (£(s— )nC)\{x}, otherwise we set^(s)nC = £,(s— )nC. 
In the first case we say that x is erased and that the domain [x, x'] has incorporated 

the domain on its left. Similarly, if s G [0, t] equals T,„ ' with A = [x, x'] G C and 

x, x' G £(s— ), then the ring at time T„ is called legal if 

aW < M^l (3.3) 



and in this case we set ^(s) n C := (£(s— ) nC)\ {x'}, otherwise we set £(s) n C = 
£,(s— ) n C. Again, in the first case we say that x' is erased and that the domain [x, x'] 
has incorporated the domain on its right. 

We point out that C n C = if C and C are distinct clusters in W t [oJ, (]. On the other 
hand, it could be C n Vt[w, C] / 0- Let x a point in the intersection and suppose for 
example that [a, x] G C while [x, 6] G Vt[oJ, C]- Then, by definition of Wt[a;, (], one easily 
derives that the Poisson processes associated to the domains [a, x] and [x, b] do not 
intersect [0, t], while at least one of the Poisson processes associated to the domain on 
the left of [a, x] intersects [0, t]. In particular, x G £(s) n C for all s G [0, t], in agreement 
with (I3.1D . The same conclusion is reached if [a, x] G Vt[u>, (] and [x, 6] G C. This allows 
to conclude that the definition of the path {£(s)} s >o up to time t is well posed. We point 
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out that this definition is t-dependent. The reader can easily check that, increasing t, 
the resulting paths coincide on the intersection of their time domains. Joining these 
paths together we get {C(s)} s >o- 

At this point it is simple to check that, given a configuration £ G jV(d m m), the law of 
the corresponding path {£(s)} s >o is that of the one-epoch coalescence process defined 
in Section[2]with initial condition £. The advantage of the above construction is that all 
one-epoch coalescence processes, obtained by varying the initial configuration, can be 
realized on the the same probability space. Given a probability measure Q on j\f(d m i n ), 
the one-epoch coalescence process with initial distribution Q can be realized by the 
random path {£'(s, -)}s>o defined on the product space £1 x j\f(d m i a ) endowed with the 
probability measure P x Q. 

3.2. Proof of Theorem 12. 131 (i) — (ii) — (in). Before presenting the proof of Theorem 
12.131 (i) — (ii) — (Hi) we state and prove a key Lemma. 

Lemma 3.1 (Separation effect). For any x G M, any configuration Q G J\f(d m [ n ) with 
x G £ any event A in the a-algebra generated by {£(s) n (— oo, x)} s ^ & any event B in 
the a-algebra generated by (£(s) n (x, oo)} s -g t it holds 

F c (AnBn{xe m}) = Pcn(-oo,x] {An{xe e(t)})Pcn[x,oo)(B n{xe £(£)}) . (3.4) 

Proof. We set Q := ( n (— oo,x], Cr := C H [x,oo), /c := Af(x,(), J '■= N(x,Q) and 
w := AA(x, Cr)- The desired result (I3.4D is implied by the following facts (i) and (ii): 
(i) For any we(] such x G £*»(£, w) the following holds. At each time s G [0, t] one has 

£, Ce (s,ij) = ^(s,^) n (-oo,x] 

if Cj satisfies for any i < k and m G N 

T« (u) = T {i+j ~ k) (Co) , U!t j ~ k) (w) = ^ +i " fe) (w) , (3.5) 

and the same identities with T and Um replaced by T and Um ■ Similarly, at each time 

s G [0, t] one has 

^ r (s,U)) =^(S,L0) fl [X, OO) 

if Cj satisfies for any i ^ k and m G N 

T«H = T (l+ "- fc) (w) , Ug\oj) = UJt u ~ k) (Co) , (3.6) 

and the same identities with T and Um replaced by T and U m ■ 

(ii) Take Co, Co G fi such that x G £&(£,o>) and x G ^ r (t,w). At each time s G [0,£] it 
holds 

^(s,co)=^(s,Co)[J^( S ,Co) 

if w G ft satisfies (I3.5D and the same identities with T and U m replaced by T and U^ 
for any i < k and m G N, and oo satisfies (13.6D and the same identities with 7" and U m 
replaced by T and U m for any i ^ A; and m G N. □ 

We first prove the renewal property for the OCP with initial distribution Q = Ren(i/, //). 
We take the special realization of the process defined by means of the universal cou- 
pling at the end of the previous section. We concentrate on the joint distribution of the 
random variables xo(t),d\(t),d2(t), proving that they are independent and giving an 
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expression of their marginal distributions. We recall that xo(t) is the leftmost point of 
£(£), while dk(t) is the length of the k-th. domain to the right of xo(t) in £(£) . 

While di(t),d2(t) are nonnegative random variables and their Laplace transforms 
are always finite, xo(t) is a real random variable and its Laplace transform could di- 
verge. Hence, it is convenient to work with characteristic functions instead of Laplace 
transforms. Given imaginary number so,s±,S2 G iR, we have 

£ ( e -sox {t)~s 1 di{t)-s 2 d 2 {t)\ 

J2 E Q {e- sox °V- s ^- s ^;x (t) = x io (0) ; Xl (t) = x h (0) ; x 2 (t) = x i2 (0)) 
io<n<*2GN 



£ j Q(dOe- SoX >o-^n-^ )-^ 2 -*n)f i0tli . 2 ( C ) , 



10<H<12& 

(3.7) 

where ( = {xk : k € N} and the function /^^^(C) is defined as the P-probability of 
the event U in given by the elements lj satisfying the following properties: 

(Pi) cHt^)n(-oo,x io ] = {x io }, 

(P2) £ c (t,uj)n [x io ,x ix ] = {x io ,x h }, 
(P3) £ C (t,w) n [x h ,x i2 ] = {x h ,Xi 2 }. 

Let us now set 

Co = Cn (-oo,xj ] , Co,i = Cn \x i(V x h ) , Ci,2 = Cn [xi i; xj 2 ] , C2 = Cn [x i3 ,oo). 

Then, by the separation effect described in Lemma [37T1 one has 

4 

A>,ii,«(C) = p(w(0) = II p ( w G ° : w fuimis ( pi ')) ' 

1=1 
where 



(PI') 


e Co (t,u;) = {x JO }, 


(P2') 


C Co,1 (*> w ) = {^io^ii}) 


(P3') 


£ Cl ' 2 (t,w) = {x^XjJ, 


(P4') 


x i2 ee C2 (t,^)- 



We stress that the last factors are (-dependent, although we have omitted ( from 
the notation. In particular, the probability P(uj e Q. : a; fulfills (Pi')) depends on £ only 
through the first point xq and the domain lengths d\, d 2 , . . . , di if i = 1, the domain 
lengths dj 0+ i, . . . , d^ if i = 2, the domain lengths c^+i, . . . , di 2 if i = 3 and the domain 
lengths di 2+ i,di 2+ 2, ... if i = 4. Thinking of ( as a random configuration sampled by 
Q, all the above domain lengths are i.i.d. with law /j, and are independent from xq 
which has law v. In particular, the random variables ( — s- P(oj e ft, : oj fulfills (Pi')) are 
independent for i = 1, . . . , 4. Using the consequent factorization and integrating over 



tQ<ll<l2< 

X 
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C in (13.7D , we conclude that 

^ ( e SOXo(t)-Sldl(t)-S2d2(t)\ 

Y^ f Q(d()P(u £l!:w fulfills (P4')) 

/ Q(dC)e- so ^oP(w e!l:w fulfills (Pi')) (3.8) 

x / Q{dQe- Sl{x ^- x ^P(uj en-.oj fulfills (Pi')) 

x f Q{dC)e- S2 ^- x ^P(u e n : w fulfills (P2')) . 

By simple computations and using that Q = Ren(z^, /i), from the above identity we 
derive that 

Eq^-" * ®-' 1 * 1 ®-" 1 *®) = L t (s Q )G t (si)G t (s 2 ) (3.9) 

where 

U(s) = PRen( 5 „, M) (o g e(*))A)W E E ®n/*( e_M!n(0) ; £(*) = K(o)» > tf.io) 

n ^ 

G t (a) = E E^(e— (°) ; e(t) = {0,x n (0)}) . (3.11) 

n> 1 

Above Lq(s) denotes the characteristic function of v, while ® nj u denotes the law of the 
SPP given by n + 1 points = xo < x\ < ■ ■ ■ < x n such that the random variables 
di = Xi — Xi-i, 1 ^ i ^ n, are i.i.d. with common law \x. 

Note that in the derivation of (I3.10D one has to keep the contribution of both the 
first and the second expectation in the r.h.s. of (I3.8D . 

By similar arguments one obtains 

k 
EQ ( e -(soxo(t)+sid 1 (t)+.:+s k d k (t))j =L t ( ao )JjG t (a i ) V/fc^O, (3.12) 

4 = 1 

with the convention that the last product over k is equal to 1 if k = 0. The above 
formula implies that the random variables xo(t) , di(t) , d2{t) , . . . are all independent, 
xo(t) has characteristic function L t and c4(£) has characteristic function Gt for each 
fc ^ 1. Note that the above arguments remain valid for sq,si,... , s& ^ (and one 
speaks of Laplace transforms instead of characteristic functions), but if E(e _sox '°( )) = 
oo we get the trivial identities oo = oo. 

(ii)-(iii) We consider now the case Q = Ren(/i). Points are now labelled in increasing 
order with the convention that xq denotes the largest nonpositive point. Similarly to 
the above proof, one can show that the random variables dk(t), k / 1, are i.i.d. and are 
independent from the random variable X\(t) — xo(t). Moreover, their common law has 
Laplace transform (I3.11D . On the other hand, due to the definition of the dynamics, £ f 
must be a stationary SPP. As a byproduct, we conclude that the law of £ t is Ren (//*), 
IM being a probability measure on (0, oo) with Laplace transform (13.1 ID . The case 
Q = Reng(/i) can be treated analogously 
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It is convenient to isolate a technical fact derived in the above proof, which will be 
the starting point in the proof of Theorem 12 .14t 

Lemma 3.2. Recall that G t (s) = / [inf (^ ))0O ) e~ sx ^t{x) and L t (s) = f A e~ sx vt{x) G 
R + U iR). Then 

Lt(s) = f Re n(8 ^(0 G m)Lo(s) E E^( e —(°) ; £(i) = K(0)}) , (3.13) 

n ^ 

G t (a) = ^ E^(e— (°) ; f(t) = {0,x„(0)}) , (3.14) 



n>\ 



where ® n n denotes the law of the SPP given by n + 1 points = xq < x\ < ■ ■ ■ < x n such 
that the random variables d\ = x\ — x^\, 1 < i < n, are Ltd. with common law fj,. 

3.3. Proof of Theorem 12. 131 (iv). Suppose that Q = Ren(<5 , /«)• Then we can write 

E Q (e-^W ; G £(*)) = E Eq^ - *^ ; G £(i) ;*i(t) = x,(0)) = 

■;gn + 

E / Q(rfC)e" SXl iPc (e(t) n [0, x,] = {0, si}) > C3.15) 

where £ = {x^ : k ^ 0}. By the separation effect described in Lemma I3TT1 we can 
write the last probability inside the integrand in (I3.15D as 

P<n[o,*,](£(<) = {0,Xi})¥ (n[XitOo) (xi G C(t)) ■ 

We observe that the last two factors, as functions of (, are Q-independent. Moreover, 
for all i G N+, it holds 



QK)P C n[^oc)(^ e £(*)) = Pe(° e £(*)) • 

Therefore, coming back to (I3.15D . using the renewal property of Q and (13.14D . we get 
E e (e-*C*) |o g m) = E /Q( d OPcn[o,x i ]U(t) = {0,a;i}) = G t (s) . 

J6N+ 

By similar arguments one gets 

k 

E c (e-^=i'^(*>|0 G £(*)) = [J G *^') ' si, . . . ,s fc G R+ UiR, 

i=i 

thus concluding the proof of Theorem 12.131 (ii). 

3.4. Proof of Theorem \2A3\(v). From (l3"TT3l) and (13341) we get that L t (s) and G t (s) 
converge to L OQ (s) and Goo(s) as i — >■ oo. This implies the weak convergence to v t and 
IM to z^oo and /i^. 



ID HIERARCHICAL COALESCENCE PROCESSES 21 



4. Recursive identities in the OCP: Proof of Theorem 12. 141 

The proof is based on the identities (I3.13D and (I3.14D in Lemma 13.21 We first point 
out a blocking phenomenon in the dynamics that will be frequently used in what fol- 
lows. Due to assumption (Al'), a separation point x between two inactive domains 
cannot be erased. As simple consequence, we obtain that the points between two near- 
est neighbor inactive domains cannot all be erased: if there exists s ^ s.t. [a, b] and 
[c,d] are inactive domains (including the cases a = — oo, d = oo) with b ^ c, then 
£(oo) n [b, c] / 0. Indeed the set [b, c] n £ is non empty (since b and c belongs to it) and 
if we assume that all points in this set are killed then the last one to be killed is for sure 
a separation point between two inactive domains and a contradiction arises. We will 
frequently use this fact below. 



By Lemma [3721 we can write, for s G M + U iR, 

oo 

G^s) = X> fc ( S ) , Ms) = E® fc+1 „(e-** ^ ; £(oo) = {0,x fc+1 (0)}) . (4.1) 

fc=0 

We explicitly compute A k (s). To this aim we consider the one-epoch coalescence pro- 
cess with law P® fc+1/ x- We observe that, due to the blocking phenomenon, the event 
£(oo) = {0, Xfc + i(0)} implies that (i) k ^ 1 and the k + 1 initial domains are all ac- 
tive, or (ii) k ^ and initially there are k active domains and one inactive domain. 
Therefore, given k ^ and 1 < j < k + 1, we introduce the following events: 

F k = { 4(0), d 2 {o), ..., 4+i(o) eA}n {£M = {o, x k+1 (o)}} , 
E Kj = { di(o) g AW e {i, . . . , k + 1} \ {j}} n {^(o) .a} n {£M = {o, x fc+1 (o)}} . 

By the above discussion, it holds 

fc+i 
A k (s) = E 8fe+lM (e- s ^+^°) ; F k )l k > , +£>® fc+lM (e- s:Efe+l(0) ; ^j) ■ (4.2) 

i=i 

The exact computation of the two addenda in the r.h.s. is given in the following lem- 
mas: 

Lemma 4.1. For each k ^ 1, it holds 

^® fc+1 M e ,* k i) (jfe + 1) . (fc _ l)! • l4 ' dJ 

Lemma 4.2. For each k ^ 0, it froZds 

E^(^ +l(0) J ^) = / MfrfaOe— 1,*A [/ ^J 6 "" 1 *^ " . (4.4) 
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We postpone the proof of the lemmas in order to end the proof of Point (i) of Theo- 
rem [2231 Due to (14. ID , (I4.2D , Lemma [4J] and Lemma |4~2"1 we obtain 

^oolSj-2^ ft _|- 11 . ffc - "H 



A=l 



(k + i)-(k-iy. 



oo „ 

+ J2 fi{dx)e- sx t x ^ A 



[Jfi(dx)i 



\k 



k=0 






-E (riw^ + E(^)-^))^- (4.5) 

oo 

= -HoW - Et(— T)T " JT]7^)T]^( S ) j + GoWe^oW 
= - V ^^ + G (s)e H ^ = 1 - e*°W + G (s)e H ^ . 

3=1 J ' 

This concludes the proof of (I2.8D (and hence of Point (i) of Theorem I2.14D . Now we 



give the proof of Lemma I4TT1 and Lemma [4721 



Proof of Lemma WM From now on we work with the one-epoch coalescence process 
whose initial distribution is given by Qk+iP- 

Let us suppose that d\(0), £^(0), . . . , dk+i(0) G A: we want to understand how the 
event F^ takes place, i.e. how points xi(0), . . . ,x/c(0) are erased while xo(0) = and 
Xfc + i(0) survive. The event F^ must be realized as follows: 

(i) the first erased point must be of the form Xj(0) with 1 < i ^ k, 
(ii) after the disappearance of Xi(0), restricting the observation on the left of Xj(0), 

one sees that Xj_i(0), Xj_2(0), . . . , xi(0) disappear one after the other, from the 

rightmost point to the leftmost point, 
(iii) after the disappearance of Xi(0), restricting the observation on the right of 

Xi(0), one sees that xj + i(0),x; + 2(0), . . . , Xfc(0) disappear one after the other, 

from the leftmost point to the rightmost point, 
(ii) and (iii) follow from the blocking phenomenon and the fact that the disappearance 
of Xi(0) creates an inactive domain, [xj_i(0),xj + i(0)]. Since the initial configuration 
has a finite number of points, the coalescence process can be realized as follows: each 
domain of initial length d waits independently from the other domains an exponential 
time of parameter X(d), afterwards if both the its extremes are still present we say that 
the ring is effective and with probability X r (d)/X(d) its left extreme is erased otherwise 
the right extreme is erased, and after this jump the dynamics start afresh. We can 
therefore describe the jumps in the coalescence process (disregarding the jump times) 
by a string a = (<ti,<72, . . . ,cr m ), where each entry <j{ is a couple <ji = (Ni,Li) with 
Ni€ {l,2,...,k + 1}, Ni ^ Nj for i / j and L, e {£, r} (N stands for "number" and L 
stands for "letter"). The meaning of ctj is the following: the domain which rings at the 
i — th effective ring is given by [xi\r i -i(0),XAr i (0)], while after its ring the erased extreme 
is the left one if L« = £ or the right one if Li = r. See figure [3] for an example. We say 
that the number Ni is associated to the letter Li. Given such a string a we denote by 
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B(a) the event that the jumps of the coalescence process are indeed described by the 
string a in the sense specified above. 



• • 



-+- 5 t/l ring: N 5 = 6,L 5 = £ 

-•- 4 th ring: iV 4 = 5, L 4 = £ 

-*- 3 rd vmg: N 3 = l,L 3 = r 

-4- 2 nd ring: N 2 = 4, L 2 = £ 

-•- l st ring: Ni=3,Li=£ 

-•- £ = 



Xq X\ 



■r-2 



x 3 x 4 



X5 



Xf, 



Figure 3. Example of a trajectory in F^, with k = 5. 



Due to our previous considerations it holds 



F k = U a 



admissible 



B(a) 



where a string a = (<7i, 02, . . . , <r m ) is called admissible if the following properties are 
satisfied: 

(PI) if L\ = £, then Ni G [2, fc + 1]; the numbers iV, associated to the letter £ are all 
the integers in [N±,k + 1] and they appear in the string in increasing order; the 
numbers iV, associated to the letter r are all the positive integers in [l,N\ — 2] 
and they appear in the string in decreasing order; 

(P2) if L\ = r, then N\ G [1, fe], the numbers iVj associated to the letter £ are all 
the integers in [JVi + 2, fe + 1] and they appear in the string in increasing order; 
the numbers Ni associated to the letter r are all the integers in [1, N{\ and they 
appear in the string in decreasing order. 

Observe that an admissible string must have k entries, i.e. m = k, and that the knowl- 
edge of (Li)i ^i^k allows to determine univocally the numbers (JVj)i <j j ^ &■ 

Recall that \}(d) = X r (d), X*(d) = Xe(d). Writing dj(0) as di (for simplicity of nota- 
tion), if a is admissible we get 



E ® fc+ i/i [(■ 



sdi(oo) . 



B{o)] = E 0fc+lAt [F(di,d 2 , -.., 4+1, 0-) 



(4.6) 
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where 

F(d 1 ,d 2 ,...,d k+1 ,a) 



Ul e M^u..,^^,,! ^^ 11 



A(di) + --- + A(4+i) A^atJ fJ£* =i A(d^.) A^) 



(the last factor is defined as 1 if k = 1). 

Observe that the law <8>fc+i/i is exchangeable, i.e. it is left invariant by permutations 
of di,d 2 ,- ■ ■ , dk+i- This symmetry leads to the identity 

^<8>k+iv[F(,di,d 2 , ■ ■ ■ , dk+i,cr)\ = M® h+lti [G(di,d 2 , . . . ,d k +i, (£j)i <i ^ fc)] 

where 

«&,*., . . . .4+,. w, « « .) = (g 1 .-*^) A(di) + AI :. (d + 'i,4 +1 ) n ^^, 

Recall that an admissible string <r is univocally determined by its letter string (-Lj)i ^ $ ^ & 
and observe that each string in {£, r}! 1 ^ is the letter string (-Lj)i <: j <; k for some admis- 
sible a. Therefore we have 

®® k+ u*[e- adlioo) ;F k ] = Yl ^ k+1 ,[F(d 1 ,d 2 ,...,d k+1 ,a)} 

a admissible 

= 2^ ^® h+1 n[G(di,d2,...,dk+l,(Li)i^is i k)] 

L lt ...,L k e{e.,r} 

= E 9k+1 „[H(d l ,d 2 ,...,d k+1 )] (4.8) 

where 

Applying Lemma [K7T1 in the Appendix with k + 1 instead of k, m = fx, f(x) = e~ sx t xf zj\_ 
and g(x) = X(x), we end up with 

E 9k+lll [H(di,...,d k+1 )] = {k - 1 _ [J ix(dx)e- sx t xsA ] k+X . (4.9) 

This ends the proof of Lemma R~T1 

□ 

Proof of Lemma \4~2\ The proof follows the main arguments in the proof of Lemma FTT1 
hence we skip some details. As in the proof of Lemma R~T1 we work with the one-epoch 
coalescence process with law ¥^ k+lll . 

Denoting the jumps of the coalescence process (disregarding the jump times) with 
the same rule used in the proof of Lemma FTTl i.e. by means of the string a, we get that 

Ekj = U CT j_ admissib i e i3(o-) (4.10) 
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where now j -admissible means that the numbers TV, associated to the letter £ appear 
in the string a in increasing order from j + 1 to k + 1, while the numbers Ni associated 
to the letter r appear in the string a in decreasing order from j — 1 to 1. Note that in 
particular a contains j — 1 letters "r" and "k + 1 — j" letters t, and therefore a has length 



fe. 



As in the previous proof we set d r = d r (0). We then compute the expectation 



iWM[e-' dl(00) ;*(*)] 



E 6 



'fc+iA* 



* 1 *«wII{ e "**' il ^6X^ fc 



*=1 



E^Afe) 



E 



®fe+lM 



^l^uwIHe-*!*^ 



KMi) 



■} 



( | e— l xftA v(dx)) E^ fcAt [Hje-^1 



AL(^) 



l di€A k 

Er=i A(rfr) 



(4.11) 



where in the second identity we have used the exchangeability of ®fc+ijU and in the 
third identity we have simply factorized the probability measure. 

Summing over j allows to remove the constraint that a must have j — 1 letters "r" 
and "fc + 1 — j" letters £, hence 



fc+i 
J2^ k+1 ,[e- sdl{oo) ;E kd ] 



r.h.s. of (I4TTT) 

fe 



( | e— 1,^ m(<**0) E®*m [n{e"**l 



^£.4 



A(di 



e;=*ak) 



(4.12) 



Applying Point (6) of Lemma IE. II (with /(x) = e ^l^g^ and g(x) = A(x)) ends the 

D 



proof of Lemma [4721 



4.1. Proof of Theorem 12. 141 (ii). The proof of Point (ii-a) is trivial, since A r = then 

xo(t) = xo(0) for any time t ^ 0. Indeed the first point Xo(t) of £(t) cannot be erased 

from the left due the infinite domain, and from the right due to the assumption A r = 0. 

We now concentrate on point (ii-b). Due to (I3.13D . we can write (for sR + U G iR) 



Loo(s) = PRen(5o,.)(0 e e(oo))Lo(a)5^ J3 fc (s) 



(4.13) 



fc=0 



where B fc («) = E^e"*^ ) ; £(oo) = {x fc (0)». 
Lemma 4.3. B (s) = 1 while, for any k ^ 1, it holds 



Bk(s) = E 



II' 






-sc(, 



^e(di)ldieA 



(4.14) 

Proo/ We work with the one-epoch coalescence process with law P® feAt . The case k = 
is trivial. We take k ^ 1. Due to the blocking phenomenon, there is only one possible 
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way to realize the event {£(00) = Xfc(O)}: only the points xo(0),xi(0), . . . ,£fc_i(0) must 
disappear, one after the other from the left to the right. Setting di = dj(0), this implies 
that d\,...,dk belong to A. In this case, knowing £(0), the above event has probability 

AJ(dl) „ A|(d 2 ) w x \*Ad k ) _yr \}{di) 



E -=i A(d,) E -= 2 M) A (4) f = \ zU x(dj) ' 

Since x k (0) =d x + d 2 ^ h 4, we get ( I4.14D . □ 

Since A^ = 7A r we have A* = 7A|. In particular A = A| + A* = (1 +7)A|. Hence, due 
to (I4.13D and Lemma 1431 we get 

where C := PR e n(5 ,^)(0 G £(°°)) an d where, in the last series, the addendum with 
k = is defined as 1. Applying Point (b) of Lemma lETTI (with f(x) = e~ sx t xe ^ and 
g(x) = A(x)), and recalling that Hq(s) = J e~ sx t x( z^fi(dx), we end up with 

Since Lo(0) = Loo(0) = 1, the latter identity applied to s = leads to C = exp < — °^ ^ > 
which in turn leads to Q2.9D . Then (12.10D follows immediately by noticing that 
PRen(i/, M )(a;o(0) G £(00)) = PRcn(5 ,M)(° G f (°°)) and from the above definition of C. 

5. Analysis of the recursive identity (I2.8D in OCP 

As mentioned in the introduction, a crucial tool to prove Theorem 12.191 is given by 
a special integral representation of certain Laplace transforms, which makes the iden- 
tities (I2.8D and ( I2.9D finally treatable. We first consider (I2.8D , focusing our attention 
on the one-epoch coalescence process in the same setting of Section [2] (i.e. the active 
domains have length in [dm in, dmax))- In what follows, we present an overview of the 
global scheme, postponing the proofs at the end of the section. It is convenient to work 
with rescaled random variables. More precisely, in the same setting of Theorem 12.131 
we call Xq, Xqo some generic random variables with law fj,, n^, respectively. Then we 
define 

Zq = ^o/dmin and -^oo = ^oo/dmax 

as the rescaled random variables. Setting for s > 



g (s) = E(e- sZo ) , 9oo (s) = E{e~ sZ ^) , h (s) = E(e~ sZ °; Z < a) 



"max 
^min 



(5.1) 



equation (12.8D becomes equivalent to 

l-goo(as) = (l-g ( S ))e h0 ^. (5.2) 

By definition and because of assumption (A2) we have Zq > 1, Z^ > 1 and a G [1,2]. 
These bounds will reveal crucial later on. 
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For later use, we point out some simple identities. We recall the definition of the 
exponential integral function Ei(s), s > 0: 

co e -t 



Ei(s) = / — dt = / dx . 

Js t Ji x 

Given a Radon measure t on [0, oo) (i.e. a Borel nonnegative measure, giving finite 
mass to any bounded Borel set), by Fubini's theorem it is simple to check that 

/oo — s(l+x) roo />oo 

— t(dx) = / due~ u / e~ ux t(dx) . (5.3) 

Above and in what follows, we will write / 6 °° instead of J, , for c ^ 0. If t(dx) = codx, 
the quantity in (15.31) is simply the exponential integral Ei(s) and the r.h.s. of (15.31) gives 
an alternative integral representation of Ei(s). In particular, the limit points in Theorem 
|2. 191 have Laplace transform of the form 

/-co -s(l+x) . - /-oo /-oo 

g^\s) = 1 - exp I - / — t{dx)j = 1 - exp | - / due~ u / e~ ux t{dx)^ 

(5.4) 
where t(dx) = c$dx. 

This observation suggests to write the Laplace transforms go, <?oo in the form Q5.4D 
for suitable Radon measures to an d too- The following result guarantees that such an 
integral representation exists. 

Lemma 5.1. Let Z be a random variable such that Z ^ 1 and define g(s) = E[e~ sZ ], 
s ^ 0. Let w : (0, oo) — s> M be the unique function such that 

g(s) = 1 -exp | - / due~ u w(u)\ , s>0, (5.5) 

w(s) = - f 9 ' { *\ , S >0. (5.6) 

Then the function w is completely monotone. In particular, there exists a unique Radon 
measure t(dx) on [0, oo) (not necessarily of finite total mass) such that 

/■oo 

w (s) = / e - sx t(dx), s>0, (5.7) 

JO 

and therefore 

/*00 o ( 1 I j* | 

<?(s) = 1 - exp { - / ^— t(dx)\, 0. (5.8) 

*- Jo i- + x > 

Moreover, 

limsup- / 9 V. g [0,11. (5.9) 

We recall that a function / : (0, oo) — >• E is called completely monotone if it possesses 
derivatives D n f of all orders and 

(-l) n D n f(x) ^0, Vx>0. 



i.e. 
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Due to the above lemma, there exist and are univocally determined two Radon mea- 
sures to and ^ on [0, oo), such that go and g^ admit the integral representation (I5.8D 
with t replaced by to and too, respectively. 

In order to rewrite (I5.2D as identity in terms of to and too, we need to express the 
function h$ in terms of to- The following result gives us the solution: 

Lemma 5.2. Let Z be a random variable such that Z ^ 1 and let g(s) be its Laplace 
transform. Let t be the unique Radon measure on [0, oo) satisfying (15.8ft and call m(dx) 
the Radon measure with support in [1, oo) such that 

m (A) = / ^^t(dx) . (5.10) 

Jo 1 + x 

For each k ^ 1, consider the convolution measure mS k > with support in [k, oo) defined as : 

/OO /"OO /"OO 

m(dx 1 ) I m(dx 2 )--- m(dx k )l Xl+X2+ ,„ +Xke A ■ (5.11) 

Then the law of Z is given by 

00 i -nfc+i 



k=\ 
In particular 



, -s(l+x) 



E\e~ sZ ;Z <a\ = t(dx) , OO. (5.13) 

J[0,a-1) 1 + x 

We point out that, given a bounded Borel set A, the series m*(A) = YlkLi jp — m^(A) 
is a finite sum, since m^ has support in [k, oo). The thesis includes that this sum is a 
nonnegative number and that the set-function A h-> m*(A), defined on bounded Borel 
sets, extend univocally to a Radon measure on all Borel sets. 

The above equation (15.13ft allows to write ho(s) in terms of to- Collecting the above 
observations we get for s ^ 

9o(s) = l-exp|- / 1 t (dx)j , 

5oo(s) = 1 - exp <^ - / — — too(dx) \ , 

I J 1 + x ) 

-s(l+x) 



h (s) = / — — t (dx) . 

J[0,a-1) l + x 

Due to the above identities, (15.21) is equivalent to 

oo g— as(l+x) i- „— s(l+x) 

t OQ {dx)= I — ■ t (dx), s^O. (5.14) 



1 + X J[a-l,oo) 1 + X 

It is convenient now to introduce the following notation. Given an increasing func- 
tion 4> '■ [0, oo) — > [0, oo) and a Radon measure m on [0, oo), we denote by m o <f> the new 
Radon measure on [0, oo) defined by 

mo(f>(A)=m((j)(A)), AcRBorel. (5.15) 
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Note that m o cj) is indeed a measure, due to the injectivity of <f>. Moreover, it holds 

/*oo r 

/ f(x)mo ( j ) (dx)= f(<THz)M<te)- (5.16) 

JO J[<f>(0),oo) 

We are finally able to give a simple characterization of (I5.14D , which we know to be 
equivalent to (I5.2D : 

Theorem 5.3. Consider the linear function <f> : [0, oo) — )• [0, oo) defined as <j>(x) = a(l + 
x) — 1. Then, equation (15.14D (and therefore also (I5.2D ) is equivalent to the relation 

too = (l/a)too<f>. (5.17) 

5.1. Proof of Lemma f5. 11 First we prove that w is a completely monotone function. 
Since g(s) < 1 for s > 0, we can write w; = f^2^L g k where f(s) = —e s g'(s). Triv- 
ially, g is a completely monotone function. Since the product of completely monotone 
functions is again a completely monotone function (see Criterion 1 in Section XIII.4 of 
HFe2IP . we conclude that g k is a completely monotone function. Since the sum of com- 
pletely monotone functions is trivially completely monotone, we conclude that ^2"^L g k 
is completely monotone. It remains to prove that / is completely monotone. To this 
aim we observe that by the Leibniz rule 

D n f(s) = -£ (?WVW«) = -^jrffjD^gis) 

= -e s J2 (?) (-l) k+l ^(e~ sZ Z k+l ) = e s E(e~ sZ Z f^ (?) (~Z) k ) 

k=0 ^ ' k=0 ^ ' 

= e s E(e~ sZ 'Z{1- Z) n \ 

Since 1 — Z ^ 0, the sign of the n-th derivative D n f is (— l) n . 

At this point, we can apply Theorem la in Section XIII.4 of HFe2H to get that there 
exists a Radon measure t(dx) on [0, oo) (not necessarily of finite total mass) satisfying 
(I5.7D . Moreover, the above measure t is univocally determined due to the inversion 
formula given in Theorem 2, Section XIII.4 HFe2H . Finally, we derive (15.8ft for s > from 
(I5.3D . (15.51) and (15.7D . The extension to s = follows from the Monotone Convergence 
Theorem. 

In order to prove Q5.9D we observe that ye~ y ^ 1 — e~ y for all y ^ 0, thus implying 
that 

-sg'(s) = E(sZe~ sZ ) < 1 - E(e~ sZ ) = 1 - g(s) , Vs > . 

In particular, the ratio in Q5.9D is bounded by 1. On the other hand —sg'(s) = E(sZe~ sZ ) > 
while 1 — g(s) > 0, thus implying that the ratio in Q5.9D is positive. 



5.2. Proof of Lemma [5721 Due to the definition of m(dx), we can write 

/■oo 

-t(dx) = / e- sx m(dx) . (5.18) 



I 



oo — s(l+x) fOO 



1 + X 
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By (I5.8D , since g(s) < 1 for s > 0, we get that the above quantities are finite as s > 0. 
Using the series expansion of the exponential function we can write 

i-«p{-/ - rj _ <(< b)} = j;L_J_(/ .-,„(*)) . ( 5.i9) 



Since 



/ e- sx m(dx)J = / e- sx m {k \dx), (5.20) 

we can rewrite A5.19D as 

E^W e-mW(dx)=X:(E^i) C5.21) 

where 

o fc j = ^— y e- M mW(dx) , /, = [j,i + 1) for j > 1 . 

Using again the series expansion of the exponential function and also (I5.20D . we con- 
clude that 



oo oo oo 



J^ ^ |o A fc | = E ry / e~ sa; m (fe) (cte) = exp { / e" sa; m((ix)} - 1< oo . (5.22) 

fe=i j=i fe=i ' ^° ^° 

In particular, we can arrange arbitrarily the terms in the series given by the r.h.s. of 
(I5.21D . getting always the same limit. This fact implies that 



oo 3 oo } 



r.h.s. of (I5T2T1) = E ( E afc J lfc odd ) + E ( E ak >i lk even ) (5 - 23) 



3=1 k=X j=l fc=l 

/*oo /'OO 

= / e-* x v + {dx) - / e" sa; z/_(dx), 
Jo JO 

where the Radon measures u + and z/_ on [0, oo) are defined as follows 

* + (A)=f;^m«(A), 
fc=i 

oo 

fe=l 

We point out that for any bounded Borel subset A C [0, oo) the above series are indeed 
finite sums since each m^ has support in [k, oo). In addition, z^ + and z/_ have support 
contained in [1, oo) and [2, oo), respectively. 

Collecting (f5T8l) . (f5TT9l) . (fOD and (f533l) . we obtain that 

/■OO /'OO 

ff ( s ) = / e - sx i/+(dx) - / e" sx ^_(dx) 
Jo Jo 

for all s > 0. Writing pz for the law of Z, the above identity implies that the Laplace 
transforms of the measures pz + v~ and v + coincide on (0, oo). Due to Theorem 2 in 
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Section XIII. 4 HFe2H . this implies that p z + V- = v + . It follows, that 

p z (A) = v+(A) - v~(A) , Vicl bounded and Borel . 

Since for A as above we can write v+{A) — v-{A) = Y^k=i /fci — m^ k \A), we get that 
the law pz coincides with (I5.12D . 

It remains now to prove (I5.13D . To this aim we observe that, since mS k ^ has support 
contained in [k, oo), the measure (I5.12D equals m on [1, 2). Since a<2 and using the 
definition of the measure m given by (I5.10D . we obtain that 

p p p — sil"4~icj 

E[e- sZ ; Z < a] = e- sx p z (dx) = e- sx m(dx) = — t(dx) . 

J[l,a) J[l,a) J[0,a-1) 1 + x 

This concludes the proof of (I5.13D . 

5.3. Proof of Theorem 15.31 We write p{dx) for the measure in the r.h.s. of (I5.17D . 
Using that a[<fi~ l (x) + 1] = 1 + x, we obtain for s ^ that 

g — as(l+x) i- g— s(l+x) r g — s(l+x) 

p(dx) = a^ 1 / _ — rt (dx) = / — — t (dx) . 



'0 l + x J[<t>(0),oo) a 1 ( 1 + a; ) " 7[a-l,oo) 1 + 3! 

The above identity implies that (I5.14D holds if and only if 

/■oo —as(l+x) /-oo — as(l+x) 

/ — t 00 (dx)= — — p(dx), VO0. (5.24) 

Jo 1 + x Jo 1 + x 

We write moo and m! for the measures on [1, oo) such that 

rriooiA) = l 1 1+ , xe too(dx) , m'{A) = \ 1+3:g p(dx) 
Jo 1 + x J l + x 

for bounded Borel subsets A C [l,oo). Then, by (I5.24D , we get that (I5.14D holds if 
and only if the Laplace transforms of the measures m^o and m' coincide on (0, oo). 
By Theorem 2 in Section XIII. 4 IIFe2ll . this last property is equivalent to the identity 
moo = to', which is equivalent to too = p- 

6. Hierarchical Coalescence Process: proofs 

6.1. Application of the recursive identity (I2.8D to the HCP. 

We begin by collecting some useful formula for the hierarchical coalescence process that 
we derive from results obtained for the one-epoch coalescence process in the previous 
section. These formula will be used throughout the whole section. 

We use notation and definitions of Theorem 12.191 In particular p and v are prob- 
ability measures on [l,oo) and R respectively. We define here X^ n \ n e N + , as the 
length of the leftmost domain inside (0, oo) at the beginning of the n-th epoch, i.e. 
x (n) = x ^\o) _ x j n )(o). Moreover we set Z^ = X^/d^. Note that X^ has law 
p( n \ Also, E stands for the expectation with respect to the hierarchical coalescent pro- 
cess starting indifferently from Q = Ren(v,p), Q = Ren(p) or Q = Ren^(p). For any 
n G N + and any s ^ let 



g (n) 



(s) = E (e"^ M ) , h^(s) = E (e-'^l, < zM<an ) (6.1) 
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where a n = df- n+l ^ /S- n >. Thanks to Theorem 12. 13[ (see also (15.2D ). we get a system of 
recursive identities 



1 -j^K-O = (1 - 5 „_ 1 ( S ))e h(n " 1) W Vn ^2 



(6.2) 



These recursive identities will be essential in the subsequent computations. Since 
Z( n ) ^ 1, by Lemma I5TT1 there exists a unique measure t^ on [0, oo) such that 



g (n) (s) 



1 — exp 



oo -s(l+cc) 



1 + X 



-t 



n) (dx)}, 



n>\. 



Invoking now Theorem 15 .31 we conclude that 



n > 2. 



(6.3) 



(6.4) 



where (j) n (x) = a n (l + x) — 1. 

Up to now we have only moved from the system of recursive identities (I6.2D to the 
new system (I6.4D . But while the former is highly non linear and complex, the latter is 
solvable. Indeed if we define 



1p n (x) := (pi o 2 o • • • o (j) n (x) 

then ip n (x) = S n+1 \l + x) - 1 and (15.15ft together with (f6~4l) imply 

1 



f 



(n) 



dW 



t (1) O Vn-1 



n ^ 2, 



Finally, using (15.131) and (16.61) . it is simple to check that 



hF>{i 



,(14*)/dM( 1 + a .)-l t (l) (da .) 



n^ 1 



(6.5) 



(6.6) 



(6.7) 



'[dO)-l,d("+ 1 )~l) 

where we used the identity (1 + ^"^(x)) = (1 + x)/d^ n \ 

6.2. Asymptotic of the interval law in the HCP: proof of Theorem 12.191 

S ection IBTTI provides us with most of the tools necessary for the proof of Theorem 12.191 
In particular our starting point is the identity (16.31) : 



Defining 



5 (™)(s) = 1-exp 



U {n) (x) 



OD -s(l+x) 



1 + X 



-t 



(n) 



(dx)}, 



n>\. 



t (n) ([0,x]) ifx^O, 
otherwise . 



(6.8) 



(6.9) 



we get that U^ is a cadlag function, dU^ = t^ and U^ n \x) = for x < 0. By (IBTBI) 
it holds that 



U in) (x) 



1 



1 

dH L 



tf (1) ty„-l(z))-tf ( %„-l(0)-)' 

[/(l)^W(l + x ) _ 1) _ [/(I) ((rf(n) _ 1)_) 



(6.10) 



n>\. 
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If we fix n > 2, integrate by parts and use 1/^(0— ) = 0, we can rewrite the integral in 



6M as 



oo -s(l+x) roo -s(l+x) 



1 + X In 1 + X 



1 "f x Jo 



tW(dx) = / — dU {n) {i 



e -s(i+y) f°° d e ~ s{ - 1+x ^ 



= i im 1 U^ly) - -f- C/^(a;)(fx. 

(6.11) 

We now use the key hypothesis (I2.14D . Since g^(s) = g(s) because d^ = 1, if it/ 1 ) 

denotes the Laplace transform of t^ (i.e. w^(s) = J °° e~ sx t^\dx)), then (I2.14D 

together with (15.61) implies 

\[msw il) (s) = c . (6.12) 

40 

Finally Tauberian Theorem 2 in Section XIII. 5 of HFe2H shows that (16.121) gives 

lim ^ = c . (6.13) 

y^oo y 

The above limit together with (16.101) implies that there exists a suitable constant C > 
such that 

U {n) {x) <C(l + x), n>l,ar>0. (6.14) 

In particular, the limit in the r.h.s. of (16.1 ID is zero and 

roo -s(i+x) P°° ( d ( e~ s ( 1+x ^\\ , \ 

I — ''"'^--l (^(TT^)) C " ,,W<fc ' "> 2 - t6 ' 15) 

By (16.101) . (16.131) and the fact that c n ->• 00, lim™-^ £/( n )(x) -> c x for all x ^ 0. This 
limit together with (I6.14D allows us to apply the Dominated Convergence Theorem, to 
get that 

roo e -s(l+x) f°° ( d ( e~ S ( 1+X ^>X\ roo e -s(l+x) 

lim / t^ n >(dx) = -c / -^ xdx = c dx 

n^ooj 1 + X V ' U 7o \dx\ 1 + X J J U 7o 1 + X 

(in the last identity we have simply integrated by parts) . In conclusion we have shown 
that g^ converges point-wise to the function g^' defined as in (12.15D . Since in ad- 
dition lim^o gco ( s ) = 1j by Theorem 2 in Section XIII. 1 of [|Fe2ll . we conclude that 
is the Laplace transform of some nonnegative random variable Zcq and that Z^ 



(00) 



9c a 

weakly converges to Z^' . 



r (oo) 



Finally, Lemma 15.21 allows to determine the law of Zc . Indeed, the measure 
m associated to t(dx) := co^x by means of (15.10D is simply of the form m(dx) = 
(cq/x)1 x j> ±dx. In particular m^ k \dx) = CQPk{x)t x j> ^dx with pj, defined in (12.161) . It 
remains then to apply (15.12D . 

Remark 6.1. It is useful to observe that if the initial scale d^ 1 ' was different from one than 
necessarily g(s) ^ g^(s). However, and that is the reason why we could fix d^ = 1, 
the limit (12.14D is invariant under rescaling the variable s by a constant i.e. (12.14D for g 
implies the same limit for g^\ 
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6.3. Asymptotic of the first point law in the HCP: Proof of Theorem 1 2. 241 We first 
prove the result for the special case v = 8q. We set 



l {n \t 



E 



exp 



{-sX^/d^} 



s <E 



Recall the notation of Section I6TT1 and in particular the definition of the constants a n = 
^(n+i)^(n) gy applying to each epoch the key identity ( I2.9D we get the recursive 
system: 



£ (n) {s) = ^-D (s/an _ l) exp {^_[/ l (n-D (s/an _ l) _ ^(n-1) 

1 + 7 



n > 2. 



Since OjOj+i • • • a n -\ = S n > /$i> by combining the above recursive identities we get 



£(n) ( 



6 



n-1 

+ 7 3=1 

We now use the integral representation (16.7ft to get 

h^\sS j) /S n) )= f {l + x)- 1 e"^ (1+x) t^(dx), j>l. (6.17) 

J[dU)-i,dU+i)-i) 

This allows to write 

n—l ., 

Setting U^ (x) = t*- 1 -* ([0, x\), we can use integration by parts and the change of variable 
y = (1 + x)/d^ to conclude that 



(l + x)- l e- s ( 1+x V din) tU(dx). (6.18) 






dW 



-s.V 



[l/<i (n) ,l) 



•s 




1 1 


— 


+ 


y 2 


IV 





uW{dWy-l) 



dy. 



In particular, we can write 



F<»)( s )-F<»)(0)=(^-1)- C,<1,(d< " l " 1) 



+ 



•ST 



-sy 



d(n) 
^W(d(»)y-1) 



dy + 



-iv-lU^id^y-l) 



dy. 



[i/d(«),i) rf(n) y i[i/d(»),i) y rf (n) y 

We have already observed that (I2.14D together with a Tauberian theorem implies the 
limit (I6.13D . Since d^ — > oo, we can then apply the Dominated Convergence Theorem 
to conclude that 

lim fWfi) - fW(0) = c fe~ s - 1 + / se-^dj/ + / 6 ^ ~ l dy] . (6.19) 
n ^°° v 7(0,1) J(o,i) V ' 

Collecting (16.161) . (16.181) and Q6.19D , we conclude that for any s & M+ the sequence 
(.£W(s)) > ± converges to 



exp 



co 



1 + 7 7(o,i) 



1 - e~ sy 



-dy 
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Since the latter is continuous at s = we get the desired weak convergence (cf. Theo- 
rem 3.3.6 in flU). 

Now we prove the result for a general v. By translation invariance, for any x G R, 

Xq = x ) = IPRen(5 u) ( Xq = ) . Hence, for any bounded continuous 



function /, 



E- 



Ren(M (f{X™ /<&>)) -- j u(dx)E^ n{5x ^ (f(X^/d^)) 

z,(dx)E Ren(5oiAt) (/((xi n) -x)/ciW)) 

and the result follows from the case v = 5q considerer above once we use the assump- 
tion linin^oo d^ = +00. This concludes the proof. 

6.4. Asymptotic of the survival probability: Proof of Theorem l2.25l This section is 
dedicated to the proof of Theorem l2.25[ We use the notation and definitions of Section 
IBTfl We start with Point (i). 

6.4.1. Proof of (i). As in the proof of Theorem 12.241 it is enough to consider the case 
v = 5q. Recall the definition of fi^ introduced before Theorem 12.191 By a simple 
induction argument based on Theorem 12. 131 ("iD. if the initial law Q is Ren(<5o, /u) then 
the law of £^')(0) i.e. the SPP at the beginning of the j th -epoch, conditional to the event 
{0 G C (0)} i s Ren(5o ,1^ )• Hence, by conditioning and by using the Markov property, 
we get 

F Q (X^ +1) = o) = P S (X« = 0) f[ P S (X? +1) = | X® = 

i=i 

= Yl F Ren(5 ,^)){ X +1 = °) ■ 

In the last line, we also used the trivial equality Pg [Xq = 0) = 1. Theorem 12. 141 ( 



ensures that 



in 



ft(j)(0) 



Wo,A««>) W +1) = 0) = e"— , Vj > 1 



where, thanks to (I6.7D , 

fc<*)(0) = ^([dW^Ci+i))) = f (l + x)-H {1) {dx) . 

J[dU)-l,d(i+ 1 )-l) 

It follows that 

F S W" +1) = 0) = exp { - -J- V ^)(0)) = exp { - -i- / ^). 

(6.20) 

If [/"^(x) = i^QO, x\) and using integration by parts one gets 



L 



t^(dx) _ C/( 1 )(d( n+1 ) - 1) /• U<U(x) 



[o,d(™+ 1 )-i) 1 + s d(™ +1 ) 7 [0)d (n + i)_ 1) (1 + x) 
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As in (16.131) our assumption implies that lim^oo ^ = cq. Since S n ^> — > oo we get 

immediately that uW ^+u~^ = c + o(l). On the other hand, if A = v /ln(d( n + 1 )) and 



d(»+i) 
using again that lim^oo ^ = cq, we have 



/ 

JfO. 



' dx = K —+dx + / — — K —+ — —dx 



'[o,rf»+*)-i)(l + x) 2 J[o,A) (1 + x) 2 J[A,d(^)-i) co(l+x)l + x 

^U<V(A) + (l + o(l)) I -^dx 

J[A,d^+ 1 )-l) 1 + X 
= (1 + O (l))coln(d( n+1 )). 

Similarly 

f £%d*> [ ^^-^dx = (l + o(l))c Hd^). 

J[0,d^+ 1 )~l) [1 + x ) J[A,d(" + 1 )-l) co(l + x) I + X 

In conclusion J [od( ra+1)l) ^-^- = (1 + o(l))c ln((i( n + 1 )). Result (») of Theorem [Z25l 
follows from (16^201) . 

6.4.2. Proof of (ii). The second part of Theorem 12.251 follows from part (i) using the 
universal coupling introduced in Section [37T1 

We distinguish between two cases. Assume first that 7 = 0. This implies \i = 0. 
In turn, site cannot be erased from any ring of its left domain. Hence, the event 
{0 e £W (°°)} depends only on the rings of the domains on the right of 0. Therefore 

FW(HO) (0 G e (n) (00)) = P R en(,W) (4" +1) = 0' 

and the expected result follows at once from Point (i) (with 7 = 0). 
Now assume that 7 > 0. Then, by Lemma [37T1 we can write 

IWlo) (0 e e (n) (oo)) = F^ en(5()iAt) (4" +1) = 0) x P Ren(5ojAt) (0 G £ (n) (°°)) 

where P* denotes the probability measure with respect to the hierarchical coalescent 

process built with Af 1 = \[ n and AJ, n = A r (i.e. the mirror with respect to the 

origin of the hierarchical coalescence process built with Ar and Xy" '). The identity 

AJ, = 7A r n implies A^ n = -Ar * . Hence, by applying twice the result of part (i), 
we get 

= (i/rfi" +1 > 

and the proof is complete. 
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6.5. Convergence of moments in the HCP: Proof of Proposition 12.221 The proof 
of Proposition 12.221 will be divided in various steps. First we will prove the result 
for f(x) = x k , and then for a generic function / satisfying |/(»)| < c(l + x k ). The 
parameter k ^ 1 is fixed once for all. 

In what follows we will use the notation and the definitions of Theorem 12.191 In 
particular p, and v are probability measures on [1, oo) and R respectively, X^ is a ran- 
dom variable with law u^ chosen here as JfW = x^\o) - x{ n) (0), Z^ = X^/a n ^ 
and , Z(°°) = ZJ oo) is the weak limit of Z^ proven in Theorem [27T91 Recall that 
^(n) _ a n-i an( j j n p ar ticular d^ = 1. Also, E stands for the expectation with re- 
spect to the hierarchical coalescent process starting indifferently from Q = Ren(z/, //), 
Q = Ren(/i) or Q = Reng(/i). Following Section [6\2l for any n ^ 1 and any s > 

we introduce g^ n \s) = E (e~ sZ " J, the Laplace transform of Z^ n \ and h^ n \s) = 
E (e~ sZ(n) l 1 < Z (n )<a ). Thanks to Theorem [2J]3] (see also (Ex2l) ) we have 



l-<7 (n) (as) = (l-g (n - 1) (s))e M " 1)(s) Vs ^ 0, Vn ^ 2 . (6.21) 

Notation warning In the sequel for any pair of C°° functions /, g the symbol D k f(x) 
will stand for k th -derivative of / computed at the point x while the symbol D k f(g(x)) 
will denote the k th -derivative w.r.t x of f(g(x)). 

The above recursive identity will be very useful in our computations. Note that by 
Lebesgue's Theorem, for any n and any k, E([Z( n ) ] k ) = l\m s ^o(—l) k D k g( n \s) € [0,oo]. 
We shall write, for simplicity, D k g^(0) := lim^o D k g( n \s). It is not difficult to prove 
by induction on n that \D k g( n \o)\ < oo, by taking the ^-derivative of both sides of 
(16.211) . using Leibniz rule and the fact that E([Z (1) } k ) < oo. In turn 

E([z(">] fc )<oo Vn^l. (6.22) 

As a technical preliminary we prove that the above bound holds uniformly in n. 

Lemma 6.2. Assume that fi as finite k th -moment, i.e. E([Z^ ] k ) < oo. Then 

supE([Z^] k )<oo. 

n > 1 

Proof of LemmaW2\ It is not restrictive to take n > 2. By (16.22D . E([Z^] k ) is well 
defined. Moreover, E([Z^] k ) = (-l) k D k g^(0). Hence, since x k e~ x sC B := k k e~ k 
for x ^ 0, we have 

D k+1 g { - n \u) du. 
The above bound and Lemma 1631 below imply that 

du 



E([Z<»> ] k ) < U + B + - 2e ° n _i /' D*+V 1} 
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for some constant A that depends on k and on E([ZW ] k ). Now by definition of g^ 
and Fubini's Theorem, we get that 



DM ° m (jft) i * - e [sy i] i' +i ^ {-f^ 



da 



,n-i]E r^(i) I* i 



exp 






<a 



n-1 E ([Z^ 1 



)ifc 



Therefore, 



EffzW ] fc ) < -A + B + -^-E ( [ZW ]* 
2 a — 1 V 



and the expected result follows. 



□ 



Lemma 6.3. Assume that /i has finite k-th moment. Then there exists a positive constant 
A (that depends on k and E([Z^ ] k ) but does not depend on n) such that 

lea 



(a — l)a n 
Proof of Lemma [631 Iterating (16.2 ID we get 



D k+1 g^ 



u 



Vn ^ 1, Vu > 0. 



n-l 



1 _ g (n + D {s) = (i_ g W (JL) ) exp I J- ^+1) (_f_) } Vs ^ 0, Vn ^ 1. (6.23) 



j=0 



Hence, by Leibniz formula, 



fc+i 



D fe +V n+1) ( S ) = E 



£=0 



fc+1 



d-1 



D k s +l -y i) -1 ^(ex P {E^ +1) (-^)} 



i=o 



^ n— 1 



+ £ 

£=i 



fc + l\ 1 



£ / n(fe+l-i) 



j=0 



D k+X-l g {l)' 



n-l 



xWz^te 



3=0 



n-l 



+ 9 



7 (D 



|^(ex P {E^)(^-)}). (6.24) 

i=o 

In order to bound L>f (exp {E"=o ^ (i+1) („5=t) }) one has to deal with 

n — l n—X -. 

y D ih^ (—\ = y -^—D e h^ + ^ 



j=0 



3=0 



a n -i 



By definition of h^ +1 \ we have \D e h^ +1 \s)\ = E ([Z^ +l )fe- sZU+1) t x 



Therefore 

n-l 



n-l 



E* 6 1i)hE(i)' ; Ef«-^ 



j=0 



a n-3 J\ - A^ \ a n-3 

3=0 3=0 



<-l 



<Z(i+ 1 )<a 



W> 1. 



£a*. 



+ 
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In turn, for any £ = 1,2, . . . ,k + 1, since h^ +1 ' (u) < h^ +1 ' (0) for any u and any j, 



n-l 



n-l 



0?»(exp{E^ +1 '(^j)})hCexp{E^ + «(^j)} 

j=0 j=0 

n-l 

< C exp { Y, h {j+1) (0) } , (6.25) 

j=0 

for some constant C depending only on k where we used the fact that for any F smooth 
enough and any £ ^ 1 it holds 

D e x e F W = p e (dF(x), D 2 F(x), ..., D e F(x)) e F{ - x) 



where Pi is a polynomial in the variables X\,...,X^ of total degree £, whose coeffi- 
cients belong to {0, 1, ...,(£ + 1)!}. 

Now observe that, for any £ = 1, 2, . . . , k, by definition of gVJ and since Z^ > 1, 

D k+1 ' e g^(u)\ = e([Z« ]*+i-V wZ(1) ) < E([Z« } k ) Vu > 0. (6.26) 

On the other hand, since 1 — e~ x ^ x for x ^ and since Z^ ^ 1, one has 



(i) M _ ! 



<— E(Z«). 



(6.27) 



Hence, by (f6724l) . (167251) . (167261) . (167271) and using the facts that a > 1 and /i^ +1 ) (u) ^ fctt+i) (0) 
for any u and any j, we end up with 

1 1 n-l 

^v- +1) ( S )| < ^+')+^Wj)^E* (i+1) (») 

for some constant C" depending only on fe and E([Z( 1 )] fc ). The expected result of 
Lemma [6731 follows from Claim [6741 below. 

Claim 6.4. For any n ^ 1 it froZds 

j=0 

Proof of ClaimWM Fix s e [0, 1]. Since zW ^ 1 and s e [0, 1], we have 



1 -" <1> (i) = 1 - E (^* z '")> 1 



e "" > 



2a n ' 



where in the last line we used that 1 — e x ^ § for x £ [0, 1]. Hence, we deduce from 
(167231) that 



n-l 



n-l 



i , i-^M _ (^(^{EMjFj)} s, ^ exp { E ,o«» (_!_)} . 



j=0 



i=o 
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Now, by definition of h^ and since e x — 1 ^ — x for any x ^ 0, 

n— 1 n—1 n— 1 



3=0 j=0 j=0 

n—1 n—1 

> - E ^ { zU+1) ^ < ^x.) + E ^ +1) (o) 

3=0 j=0 



*-** a n -i- x ^ v ' a n - l (a-l) *-*> v y 

3=0 3=0 v ; 3=0 

We deduce that 

3=0 V J 

Optimizing over s e [0, 1] (choose s = a „_~ finally leads to 



71-1 



2-exp{£^> (0) U4£K^> ,-_?» 
a n i *—< J a" (a — 1) a — I 

3=0 V ; 

The claim follows. D 



The proof of Lemma 16731 is complete. □ 

We can now prove Proposition 12.221 for the special choice f(x) = x k . 

Proposition 6.5. Assume that \i as finite k-th moment, i.e. E([Z^ ] k ) < oo. Then, in the 
same setting of Proposition ^. 22\ z(°°) has also finite k-th moment. Moreover, 

lim E([Z^] k )=E{[Z^] k ). 

n— >oo 

Proof. Thanks to Lemma [6721 sup r-CM+ M([Z^ n ' ] k ) < oo. Fix a decreasing sequence 
of positive numbers (s m ) mg pj that converges to 0. Since x k e~ SmX is a continuous 
bounded function on R + , Theorem 12.191 implies that \\m n ^ 00 'E{\Z'^>} k e~ SmZn ) = 

E ( [Z(°°) ] fc e- SmZ(oo) ) . Hence, by the Beppo Levi's theorem 



E ([^(oo) ]fc) = limEftzMlV 8 - 2 ' 00 ^ lim lim E([zW] fc e- SmZ( " )N ) <oo 



m— >-oo n— »oo 



Hence Z^> has finite fc^-moment. 
Next, for any s > 0, we write 

E([ZW ] fe ) = E ([ZW ] fc e" sZ(n) ) + (-l) fc (D k gW(0) - D k g^{s] 
= E ([zW] k e- sZ(n) ) + (-l) fe+1 f S D k+l g^(u)du. 
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Hence, thanks to Lemma [6731 
E ([2(») ]*) _ E([#(°°) ]*) < E([zW]V szN )-E([zH]fe) + T \D k+1 g( n \u)\d 



u 



+ 



E 



0]fc e -' z(n) ) -E([Z(°°)] fc ) +As(l + s) 



lea 



^ W (^T 



r/» 



(a - IK- 1 y 

where A is a positive constant that depends on k and EQZ^ ] k ) but does not depend 
on n. Note that, by definition of g^ l \ Fubini's Theorem and then the Dominated Con- 
vergence Theorem, 



a ./o 



£,*+l^(D 



!/ 



,n- 1 



,n— 1 



du 



Ojifcg-sZW 



E 



[z (l)]fc+l e -(«zW)/a»-i du 



E([Z (1) ] fc 



1-e 



-(sZ( 1 ))/o n - 1 



^0 



E 



) jfc e sz °° j w j ieri n tends to 



By applying Theorem |2.19| E 
infinity. Therefore, 

lim E([Z^] k )-E([Z^] k ) 

n— »oo 

< E('[Z( 00 )] fc e- sZ( ° 0, ) -E([Z(°°)] fe ) +As(l + s) Vs > 0. 

The proof is completed by taking the limit as s \. 0. □ 

Proof of Proposition \2.22\ Let / be such that \f{x)\ ^ C + Cx k . For any L ^ we define 
fi{x) = f(x) if \x\ < L, and /l{x) = f(L) if |x| ^ L. Note that by Proposition 16.51 
E (/(zJ(°°))) < oo. Also, \f(x) - f L (x)\ < 2C(1 + x fc )l w >L and f L is bounded by 
construction. It follows that 

E(/(2»)) - E(/(zT(°°)))| ^ |E(/(ZW)) _ E(/ L (Z< n >))| + |E(/ L (ZW)) _ E(/l(Z(°°))) 

+ |e(/ l (z( °)))-e(/(z( 00 ))) 

< 2C7E ((1 + [ZW ] k )l z{n) > L ) + |E(/ L (ZW)) - E(/ L (Z(°°))) 
+ 2C7E((l + [Z( °)] fc )l z(oo)>L ). 

Since f L is bounded and continuous, lim |E(/ L (Z( n ))) - E(/ L (Z(°°)))| = 0. On the 

other hand, taking L among the points of continuity of the distribution function of 
Z(°°\ using that x i-> x k i x< L and x \-t 1 x <l are bounded and Proposition I6.5[ we 
have 



Jim E ((1 + [ZW } k )t z{n) > L ) = Jim {l - E (l z(B)<L ) + E ([Z<»> ] fc ) - E ([Z™ } k t z(n)<L ) } 

= 1 - E (1 Z( . )<L ) + E ([^) ] fc ) - E ([ZH ] fc l z(oo)<L , 
= E((l + [Z(°°)] fc )l z(oo)> ^ 
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Therefore, 

lim E{f{Z^))-E{f{Z^)) < 4CE ( (1 + [Z^ ] fc )l z(oo) > L 

n— >-oo \ ^ 

Now, since E ([.Z^ 00 ) ] k ) < oo and by Lebesgue's theorem, the right hand side of the 
latter tends to when L tends to infinity. This achieves the proof. □ 

Appendix A. Proof of Lemma [27T71 

We provide the proof of Lemma [2 .171 Part (i) follows immediately from Lemma [5TTT 
As far as part (ii) is concerned, if the mean of \i is finite it is trivial to check that the 
limit (I2.14D holds with cq = 1. Indeed, by the Dominated Convergence Theorem, both 
—g'(s) and (1 — g(s))/s converge to the mean as s | 0. Let us now assume that the 
mean is infinite and that for some a £ [0, 1] F(x) := (J>((x, oo)) = x~ a L(x) for some 
slowly varying function L. Notice that F(l—) = 1. Let 

Z A (s) = I (e~ y - ye-y)y- a L(y/s)dy , 

J[As,oo) 

W A (s)= [ e-y y - a L(y/s)dy. 

J[As,oo) 

Using integration by parts and the change of variables y = sx, given A > 1 we can 
write 



- S g'( s ) = - f sxe- sx dF{x) 

J[l,oo) 

= - i sxe- sx dF(x) - sxe- sx F{x)\^_ + / (e 

J[1,A) J[A,oo) 

= - i sxe- sx dF{x) + sAe- sA F{A-) + s a Z A {s) = £ + s a Z A {s) , 
J\hA) 



(Al) 
where the error term £ satisfies \£\ < sA. Similarly, we can write 



1 _ g( s ) = 1 + / e~ sx dF(x) 

J[l,oo) 

= 1 + / e~ sx dF{x) - e~ sA F(A-) + s a W A {s) = £ ' + s a W A {s) , (A.2) 

J[1,A) 

where, by taking some Taylor expansion, we get \£'\ < C{A)s for a suitable positive 
constant C(A) depending only on A. Since the mean is infinite the Monotone Conver- 
gence Theorem and De FHopital rule imply that 

lim(l — g(s))/s = lim —g'(s) = oo . (A.3) 

Comparing the above limits with (IA.1D and (IA.2D we deduce that both s a ZA(s) and 
s a WA(s) must diverge as s goes to zero. In particular, the limit (I2.14D is equivalent to 
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the limit 

lim lim ^4 = a . (A.4) 

Afoo 40 W A (s) 

As proved in HFe2H f see Section VIII.9 there), L is slowly varying at oo if and only if 
it is of the form 

L(x) = a(x) exp { [* ^-dy] (A.5) 



where e(x) — > and a(x) — > c < oo as x — > oo. In particular, given 5 > 0, for any a; 
large enough x~ s ^ L(x) ^ x s . Since in (IA.5D e(x) — > and a(rr) — > c < oo, for any 
5 > there exists A > such that c/2 < a(x) ^ 2c and |e(x)| ^ 5 for x ^ A. Thus, for 
any s < 1/A the integral representation (IA.5D implies that 

I^^^^zji/Jj^ 4 ^"^^' y ^ s - (A - 6) 

We now distinguish two cases: 

• Case q € [0, 1). Choose 5 > such that a + 5 < 1, A > 1 and s ^ 1/A. The Dominated 
Convergence Theorem together with (IA.6D implies that 

/>oo 

limZ A (s)/L(l/s) = {e~y-ye~y)y~ a dy, (A. 7) 

limW A (s)/L(l/s) = e- y y- a dy. (A.8) 

40 y 

At this point, AA.4D follows from (1A.7D , (1A.8D and a trivial calculation. 

• Case a = 1. It is convenient to write 

Z A (s) = W A (s) - T A {s) , T A (s) := f e~yL(y/s)dy . 

J[As,oo) 

Then, (IA.4D follows if we can prove that 

lim lim sup *,[ = . (A.9) 

^too 40 W A {s) 

Given 5 > we take A > 1 and s ^ 1/A assuring (IA.6D . Then we can bound 

Ta(s) < 4f™e-y(y 5 \/y~ s )dy _ C 5C 



W A (s) (l/4e)fH(y-'AyS) j^y^dy (1 - (As) s ) ' 
The above bound trivially implies GA.9D . 

Appendix B. An example of interval law not satisfying (I2.14D 

We provide here an example of a law which does not satisfy (12.14D and therefore 
does not fulfill the hypothesis under which our main Theorem |2.19| holds. Furthermore 
we have numerically analysed the set of identities (I6.6D with v- 1 ' corresponding to 
this choice for the initial distribution. The results for the corresponding function U^ n \ 
displayed in Figure @J strongly suggest that in this case the measure yS- n > does not have 
a well defined limiting behavior as n — > oo. 
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Proposition B.l. Let G be a geometric random variable with parameter p = 1 — e~ x , 
A G (0, 1). Define X = e G and g(s) = E {e~ sX ), s ^ 0. Then, lim fl _> ^=Sy does 
not exist. More precisely, for any a G [0, 1) and any n e N, set s n = e _n_a . Then 
linin^oo i" 9 ; n ( =: L a exists, and a — >■ L Q is a non-constant function. 

i — 9{ s n) 

Note that the constraint A G (0, 1) is equivalent to the fact that X has infinite mean. 

Proo/ Fix a G [0, 1) and set s n = e~ n - a . Since P(G = fc) = p{\ - pf^ 1 for A; ^ 1, we 
have F(x) = F(X ^ x) = e -^ lnx ^ +x = x - x e H^x} w h ere |y| = z + 1 _{ z } is t h e ce iling 
function of z (i.e. the smallest integer greater than or equal to z), and {z} G (0, 1] is 
the fractional part of z. Note that F(e) = 1. 

Then, using an integration by parts and the change of variables u = sx, we have 

-g'{s) = E (Xe- sX ) = - f xe~ sx dF{x) = e x ~ se + f (1 - sx) e - sx F(x)dx 

J[e,oo) J[e,oo) 

A-se I >- - --' - ( U 



e v ~ se + - / (1 - u)e~ u F (-)du 

S J[es,oo) ^ sJ 

^-se + S X-1 J (1 _ u ) u -A e -« e A{ln«-ln S } du 

J fes.oo) 



'[es,oo) 

Similarly 



1 _ g ( s ) = l + J e - sx dF{x) = 1 - e~ se + s x f u~ x e- u e x{lnu - lns} du. 

J[e,oo) J[es,oo) 

Since {In u — In s n } = {In u + a} for any n, it follows that (recall that A G (0, 1)) 

-s a'(s ) Sne 1 ' 3 ^ + s x f, Jl - u) U - x e- u e x ^ nu+a Uu 



n^ool-g( Sn ) n->oo 1 _ g-s»e + g A f u -A e -n e A{ln«+a}^. 



[es n ,oo) 

J (0co) (l - n)-u- A e-"e A 'f lnu+a >(i'u 

-• L a . 



u 



1(0 oo) u~ x e- u e x V nu+a }du 

Suppose L a to be equal to 1 — C for all a. Then, by the change of variable v = u/j3 
where (3 = e~ a , we can write 

8 f n °° t>- A+1 e-^e A { ln ^^ (B + 1) f n °° ^-^e^e-- 8 ^ 1 ™}^ 

where B = B — 1. Consider now the functions / and § on D = {z e C : |z| <1} defined 
as 

POD 

f(z)= v- x+1 e- v e- zv e x{lnv} dv, 

Jo 

POD 

g( z )= / v - x e- zv e- v e x{lnv} dv . 
Jo 

By Fubini's Theorem and the series expansion of the exponential function, one gets 

that / and g are holomorphic functions on D. Hence, the same holds for the function 

H(z) = (1 + z)f(z) - Cg(z). Due to the last identity in (IB.lft we get that H is zero 

on a subinterval of the real line. Due to a theorem of complex analysis, the zeros of 
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a non-constant holomorphic function are isolated points. As a byproduct, we get that 
H(z) = for all zeD. 

Writing the power expansion of H(z) around z = and using that H = 0, we get 



oo 



\+i+n e -v e \{\nv} dv = (c - 1) / v - x+n e- v e x{lnv} dv , Vn ^ . 
'o Jo 

Note that it must be C > 1 By iteration we get 



oo 

v- x - 



\+n e ~v e X{\nv} dv ^ ((7 _ !)n „ ^ Q (B.2) 

'0 

On the other hand the above l.h.s. is larger than T(n + 1 — A) := J* °° v~ x+n e~ v dv. 
Iterating the identity T(z + 1) = zT(z) we get that 

T(n + 1 - A) = (n - A)(n - A - 1) • • • (1 - A)r(l - A) , 

which leads to a contradiction with (IB.2D . 

D 

Appendix C. Z-stationary SPPs 

Z-stationary SPPs and stationary SPPs have many common features. In this Appen- 
dix we point out some properties of Z-stationary SPPs, whose proof (only sketched 
here) follows by suitably adapting the arguments used in the continuous case. 

Let us suppose that Q is the law of a Z-stationary SPP, nonempty a.s. We derive here 
some identities relating Q to the conditional probability measure Qq := Q(-|0 G £). 
These identities are similar to the ones relating the law of a stationary SPP to its Palm 
distribution IIDVI1 , HFKASH . Since all random sets are included in Z it is more natural to 
work with the subspaces of M defined as 

M z = {£ G M : £ C Z} , (C.l) 

Ml = {£ G A/" z : G £} , (C.2) 

Moreover, we prefer to write r x £ instead of £ — x. 

Similarly to (I2.2D we get a simple relation characterizing Q by means of Qq. 

Lemma C.l. Given a nonnegative measurable function f on J\f z it holds 
, , *ifc)-i 

/ Q(dOf(0 = S(o g / Q (dO E f(r*0 ■ (C3) 

Proof. From the Z-stationarity of Q it is simple to derive for all measurable functions 

g : Mq ->■ [0, oo) and t G N + that 

/ Qo(dOg(0 = , n( l a ts [ Q(dO E 9(t v Z) (c.4) 

J t«lU G £) J ^ n(M 

Given a measurable map « : Z x A/J — >■ [0, oo), setting #(£) = ^2 xeZ v(x, £) in the above 
identity, we get 

E / Qo(dOv(x, = X / Q(rffl E E "fo'VO- ^ 5 ) 
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Figure 4. We consider an HCP with d^ = 2 n ~ 1 (the relevant choice 
to describe East model) with initial law specified in Proposition B.l and 
parameter q := 1 — p = 0.1,0.5,0.8. In the first case the limit (2.14) 
exists with cq = 1. Instead, for q = 0.5, 0.8 as proven in Proposition B.l 
the limit (2.14) does not exists. We plot here U^ n \x)/x for x = 10 as 
a function of n. The data indicate clearly that for p = 0.1 U^(10)/10 
converges to 1 as we have proven (see Theorem 12.191 and especially the 
comment below formula (I6.15D ). Instead for the other two choices of 
the parameters C/( n )(10)/10 has an oscillating behavior which strongly 
indicates the non existence of the limit for U (jl \x), hence for g^ n \s). We 
have checked that an analogous behavior occurs for different choices of 
x. Note that if we were instead considering for the same initial dis- 
tribution but a different choice of d^ (satisfying the basic hypothesis 
d^ n ' — > oo for n — > oo) we would get the same behavior. Indeed if 
we consider for example the choice relevant for the paste-all model, 
^(n) _ n ^ t .j ien foe plot of [A n ->( 10)/ 10 would exactly be the same as 
above but with n replaced by log 2 (n) in the x — axis (and in this case 
our data would cover 2 20 epochs). 



Reasoning as in the proof of (1.2.10) in HFKASH . for any measurable function w : Z x 

Z x Af z -> [0, oo) we get 

^^/Q( ( ^Hx,|/ ) r#(i/ee)=EE/ S ^W^' T ^ 1 (^^ CC.6) 



x&Ly&L 



xez j/e 
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Combining (IC.5D with (IC.6I) where w(x,y,£) = v(x,g)l(y G (0, t])/t we get 

E / Qo(dOv(x, = * E I Q( d 0v(x, t x £)1(x € . (C.7) 

At this point we take 

v(x,£) := l(x = x (T- x £))f(T- x t;) 

(if £ = we set v(x, £) = 0). Note that «(a:,r x £) = l(x = £o(£))/(£)> tmis implying 
together with (1C.7D that 

Q(#)/(£) = E / G(#M* S ^e)l(^ e e) = Q(0 € C) E / Qo(#M*, • CC.8) 

In order to understand the last integral, take £ e Mq. Then x = xo(t_,j;£) if and only if 
^ — x < xi(£). Therefore, changing at the end x into —x, we get 

*i(0-i 

E«(».o= E /(^) = E /( r *£)> veeA/?. (c.9) 

XGZ i62: 2=0 

Combining GC.8D and qC.9D we get the thesis. □ 

Taking / = 1 in (IC.3D we deduce that xi(£) must have finite mean w.r.t. Qq. In par- 
ticular, if Qo is the law of the renewal SPP on Z containing the origin and with domain 
length p (i.e. Qq is the law of Reno(yu)), then p, must have finite mean. On the other 
hand, given p probability measure on N + with finite mean, the identity (IC.3D univo- 
cally determines the probability measure Q if Qo is defined as the law of Ren(// 1 0). 
One can then prove that the so-defined Q is the law of a Z-stationary SPP and that 
Qo = Q(-|0 € £). Finally, as in the continuous case, relation (IC.3D allows to derive 
a simple description of Z-stationary renewal SPPs similar to the one mentioned after 
Definition 12 .71 We leave the details to the interested reader. 

Appendix D. Exchangeable SPPs 

We endow with the product topology the space £1 = (0, oo) N + of sequences of positive 
numbers, and we denote by B its Borel a-field. We write a generic element of fl as 
u = (u n : n e N+). Let £ n be the a-subfield generated by the events that are invariant 
under permutations of Z fixing all points x G N + with x > n. Let £ := n^ =1 £ n be the 
exchangeable a-field. Since Q is a standard Borel set, given a probability measure Q 
on there exists a regular conditional probability associated to £, i.e. a measurable 
map pq : £1 x B — > [0, 1] satisfying the following properties: 

(i) for each A € B, pq(-, A) is a version of P(A\£), 
(ii) for Q-a.e. we!], Pq(oj, •) is a probability measure on (Q, B). 

Due to De Finetti's Theorem, if Q is an exchangeable probability measure on O, then 
for Q-a.e. u the measure pq(uj, •) is a product probability measure on f2. The inverse 
implication is trivially true, hence De Finetti's Theorem provides a characterization of 
the exchangeable probability measures on Q. 
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Suppose that Q is a left-bounded exchangeable SPP containing the origin (see Def- 
inition I2.9K • By definition, Q has support on the subspace E C M given by the con- 
figurations £ G M empty on (— oo, 0), containing the origin and given by a sequence 
of points Xfc(£), k G N, diverging to +00. We can define the measurable injective map 
^ : E ->• f2, with \£(£) = uj and w n = x n (£) - x n -i(£). We call Q the measure Q o ^ _1 . 
Trivially, Q is an exchangeable measure on $7, hence we can apply De Finetti's Theorem 
and get Q(A) = J n Q(cLj)pq(uj, A) for all A G B, where pq(uj,-) is a product probabil- 
ity measure. Since trivially pq(uj, •) has support on \I/(S), the pull-back of /9q(w, •) is 
a well defined probability measure on H corresponding to the law of Ren(6o,p UJ ). As 
byproduct, we get 

Q(A) = / Q(du)Ren(5 , p u )[A] , A c A/" measurable . (D.l) 

The above decomposition allows to extend our results stated in Section [2] to right 
exchangeable SPPs containing the origin. We give only some comments, leaving the 
details to the interested reader. Consider for example the HCP starting from Q, i.e. 
£^(0) has law Q. By applying inductively Theorem 12 .131 we get that, given n ^ 1 and 
t e [0,oo], the law of ^ n \t) conditioned to the fact that e £ (n) (i) has the integral 
representation 

Q(do;)Ren( ( 5o,[ / Ua J ]S n) ) 

for a suitable probability measure [/^]t on (0, 00). 

In particular, if each p w satisfies the limit (12.14D for some constant cq(u), we get the 

following: fixed k ^ 1, the rescaled random variable [xj? (0) — x^\(0)]/a^l n , defined 
for the HCP starting with law Q and conditioned to the event {0 G C \^)}> weakly 
converges to a random variable whose Laplace transform g(°°) is given by 



gM(8)= / Q(dw)^ (w) (a), 

where g?^ » has been defined in (I2.15D . Note that new limit laws emerge in this way. 

Let us now pass to stationary exchangeable SPPs. One can formulate De Finetti's 
Theorem also for exchangeable laws on the space $7' = (0, oo) z of two-sided sequences 
of positive numbers. At the end we get that a stationary SPP, nonempty a.s. and with 
finite intensity, is exchangeable if and only if its Palm distribution Qq satisfies 

Qo(A) = / Q(da;)Reno(^)[^], A C M measurable , (D.2) 

where (i) p w is a probability measure on (0, 00), (ii) for any A C M measurable the 
map O! 3 uj — > Reno(/i w )[«4] is measurable (thus implying that the map ui — » p^ is 
measurable) and (iii) Q is the image of the law Q of the SPP under the map A/"o°° — ► 
(0, oo) z , mapping £ in (x k (Q - x k -i(€) ■ k G Z) (recall (|27T1) L 
Using (E2D and ([2721) we conclude that 

Q(A)= Q(du)Ren(p w )[A] , A C A/" measurable . (D.3) 

in' 
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The above decomposition of Q allows to extend our limit theorems to the HCP starting 
with law Q, i.e. from a stationary exchangeable SPPs. In particular, ^ n \t) will be a 
stationary exchangeable SPP for all n ^ 1 and all t e [0, oo]. In addition, for k 7^ 1, as 
n — > 00 the rescaled random variable [a^ (0) — ^fc_i(0)]/^ n weakly converges to the 
random variable z[°°' introduced in Theorem 12.191 



Appendix E. A combinatorial lemma on exchangeable probability measures 
The next combinatorial lemma has been used in Section [3j 

Lemma E.l. Let vn^ be an exchangeable probability measure on S k , S = (0,oo), i.e. 
mfc is left invariant by any permutation of the coordinates (si, . . . , s^) £ S k . Call m the 
marginal ofmk along a coordinate (it does not depend on the coordinate). Then, for any 
bounded function f : S — > R, and any bounded function g : S — > (0, 00), it holds 
(a) 

in. ( 9{si) k y\ 9(sj) -A- A E m (/) fe 



(b) 

F ffi g(fO/(gi) ^ E rn(f) k 

m fe 11 v-fc / s ~ U\ 



Proof of Lemma Wl] We will give only the proof of Point (a) which is a bit harder. The 
proof of Point (6) follows essentially the same lines, details are left to the reader. 

Since the law mfc is left invariant by any permutations of the coordinates (si, . . . ,Sf.) £ 
S k , we have 

_1_ v^ / s(s a (i)) TT__£^wl__A f( \\ 

= k ' k U*«> + ■ ■ ■ + «<m„) ii es **,) ii /(S " (,,) J 

1 uj. / ,, x ., , y^ 9(8*0-)) TT g( s *(0) 

where Sk stands for the symmetric group of {1, . . . , k}. Hence the result will follow 
from the identity 

v , / (s - (i)) f / n ^ =k-i (E.D 
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and the product structure of m^. Now we prove (IE.1D . Divide the sum in (IE.11) depend- 
ing on the value of <r(l) and a(k): 

'•- of ED = E -r-^-r-s E E S ^'^ 



,^1 «(»l) + ■ ■ ■ + »(«») £S „% : fi E*.,' 9(^0) ) 

The thesis will follow from the fact that, for any h,ik, the last sum in the latter is equal 
to 1 . Equivalently, we need to prove that, for any n, 

Sn gr^g) 1- CE.2) 

This is done by induction. Indeed, the thesis is trivial for n = 1. Assume that (IE.2D 
holds at rank n — 1. Then, 



n / \ n n / \ 

ST^ TT 9{Sq(i)) y^ V^ T-T gtMi)J 



<r(l)=il 



y^ g(fii) y^ rr s(mo) 



<r(l)=il 

Note that, by the induction hypothesis, the second sum is equal to 1 (for any i\). Hence, 
y^ tt 9{Sq(i)) y^ g(s h ) = 

J^fiE^Ovo)) fes(«i) + *(*) + ••• + *(*.) • 

This ends the proof of (IE.2D and thus of Point (a) . As already mentioned the proof of 
Point (b) is easier (only (IE.2D has to be used), details are left to the reader. □ 
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